

NASA-CR- 166529 
19840005090 


A Reproduced Copy 

OF 


/&&>, 


Reproduced for NASA 
by the 

NASA Scientific and Technical Information Facility 



NF02379 



hPR 24 1984 

MKGLEY research center 
LIBRARV, NASA 
HAI/PTON, VIRGINIA 


FFHo 672 Aug 65 



NASA Contractor Report 166529 



A Grid-Embedding Transonic Flow 
Analysis Computer Program for 
Wing/Nacelle Configurations 

1NASA-CR- 166529) A GRID-IMEEDCING TfiANSCKIC N84-131S3 

FLOH ANALYSIS COMPUTER PROGRAM FOR 
. HIUG/ilACELLE CONFIG UR A Tit) US Final Eeport 
(Lockheed-Georgia Co. , Marietta.) 109 p Unclas 

HC A06/MF AO 1 CSCL 0 1A G3/G2 43194 


Essam H. Atta 
Joseph Vadyak 

Contract NAS2-11285 
November 1983 






NASA Contractor Report 165529 


A Grid-Embedding Transonic Flow 
Analysis Computer Program for 
Wing/Nacelle Configurations 


Essam H. Atta 
Joseph Vadyak 
Lockheed-Geargia Company 
Marietta, Georgia 


Prepared for Ames Research Center 
under Contract NAS2-11285 



Nafonai Aeronaut’csand 
Scace Acm-n-strat-cn 


Ames Research Center 



0RIG5NAL PAGE IS 
OF POOR QUALITY. 

TABLE OF CONTENTS 


SECTION PACE 

I. INTRODUCTION 2 

II. COMPUTATIONAL METHOD 3 

1 . INTRODUCTION 3 

2. COMPONENT GRID TOPOLOGY AND GENERATION 3 

2.1 WING COMPONENT GRID GENERATION 3 

2.2 NACELLE COMPONENT GRID GENERATION 5 

3. GOVERNING EQUATIONS FOR THE INVISCID FLOW 8 

4. NUMERICAL ALGORITHM 12 

4.1 BASIC FINITE-DIFFERENCE SCHEME 12 

4.2 BOUNDARY CONDITIONS 14 

4.3 ITERATION PROCEDURE 16 

4.4 INTERGRID PROPERTY TRANSFER 16 

5. NUMERICAL RESULTS 17 

III. SUBROUTINE DESCRIPTIONS '. . . 40 

1. INTRODUCTION 40 

2. SUBROUTINE DESCRIPTIONS FOR THE GRGEN3 WING GRID 

GENERATION PROGRAM 40 

3. SUBROUTINE DESCRIPTIONS FOR THE NGRIDA NACELLE 

GRID GENERATION PROGRAM 41 

4. SUBROUTINE DESCRIPTIONS FOR THE LK GRID INTERFACE 

PROGRAM 44 

5. SUBROUTINE DESCRIPTIONS FOR THE TWN FLOW SOLUTION 

PROGRAM 44 

IV. INPUT PARAMETERS AND OUTPUT INTERPRETATION 50 

l. INTRODUCTION 50 

; 2. GRC.EN3 WING GRID GENERATION PROGRAM INPUT PARAMETERS . . 50 

2.1 NAMELIST GRID3D 50 

2.2 NAMELIST GRIDIN 52 

2.3 INPUT DATA FOR WING MULTIPLE STATION DEFINITION 

OPTION 54 

2.4 PARAMETER STATEMENT SPECIFICATION 55 

2.5 FILE USAGE 56 

3. GRGEN3 PROGRAM OUTPUT INTERPRETATION 56 

4. NGRIDA NACELLE GRID GENERATION PROGRAM INPUT 

PARAMETERS 56 

4.1 NAMELIST LIST1 57 

4.2 NAMELIST LIST2 66 

4.3 NAMELIST LIST3 68 

4.4 NAMELIST LIST4 68 

4.5 PARAMETER STATEMENT SPECIFICATION 69 

4.6 FILE USAGE 70 

5. NGRIDA PROGRAM OUTPUT INTERPRETATION 70 

6. LK GRID INTERFACE PROGRAM INPUT PARAMETERS 7 1 

6.1 NAMELIST NACEL 71 

6.2 PARAMETER STATEMENT SPECIFICATION 71 

6.3 FILE USAGE 73 

7. LK PROGRAM OUTPUT INTERPRETATION 73 



Vrr; 


8. TON FLOW SOLUTION PROGRAM INPUT PARAMETERS 74 

8.1 NAMELIST FLOWIH. . 74 

8.2 NAMELIST FLOWC 77 

8.3 PARAMETER STATEMENT SPECIFICATION 82 

8.4 FILE USAGE 83 

9. TON PROGRAM OUTPUT INTERPRETATION 83 

V. SAMPLE CASES AND SUGGESTIONS FOR USAGE 85 

1. INTRODUCTION 85 

2. GRGEN3 WING GRID GENERATION FROGRAM SAMPLE CASE .... 85 

3. NGRIDA NACELLE GRID GENERATION PROGRAM SAMPLE CASE. . . 85 

4. LK GRID INTERFACE PROGRAM SAMPLE CASE 91 

5. TON FLOW SOLUTION PROGRAM SAMPLE CASE 91 

6. SUGGESTIONS FOR USAGE 97 

VI. CONCLUSIONS AND RECOMMENDATIONS 103 

REFERENCES 104 


CRIGSJAL F/VvVI 
OF POOR QUALITY. 


7 


ii 


V* 



A GRID EMBEDDING TRANSONIC FLOW ANALYSIS COMPUTER PROGRAM FOR WING/NACELLE 

CONFIGURATIONS 


Essan H. Atta and Joseph Vadyak 
Advanced Flight Sciences Department 
Lockheed-Georgia Company 
Marietta, Georgia 30063 


SUMMARY 


An efficient grid-interfacing zonal algorithm has been developed for 
computing the three-dimensional transonic flow field about wing/nacelle multi- 
component configurations. The algorithm uses the full-potential formulation 
and the AF2 fully-iraplicit approximate factorization scheme. The flow field 
solution is computed using a component-adaptive grid approach in which sep- 
arate grids are employed for the individual components in the multicomponent 
configuration, where each component grid is optimized for a particular 
geometry such as the wing or nacelle. The wing and nacelle component grids 
are allowed to overlap, and flow field information is transmitted from one 
grid to another through the overlap region using trivariate interpolation. 

This report presents a discussion of the computational methods used to 
generate both the wing and nacelle component grids, the technique used to 
interface the component grids, and the method used to obtain the inviscid 
multicomponent flow field solution. Computed results and correlations with 
experiment are presented. Also presented are discussions on the organization 
of the wing grid generation (GRGEN3) and nacelle grid generation (NGRIDA) 
computer programs, the grid interface (LK) computer program, and the 
wing/nacelle flow solution (TON) computer program. Descriptions of the re- 
spective subroutines, definitions of the required input parameters, a brief 
discussion on interpretation of the output, and sample cases illustrating ap- 
plication of the analysis are provided for each of the four computer programs. 

An isolated nacelle transonic flow analysis computer program has been 
developed as a part of this investigation. The theory and usage of the iso- 
lated nacelle program is discussed in a separate report. 
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SECTION I 
INTRODUCTION 

Reliable and efficient three-dimensional transonic analysis methods are 
required to make realiotic and coot-effective predictions of aircraft aero- 
dynanics. Early efforto to predict the transonic flow field about aircraft 
multicomponent configurations are baced on the transonic snail dioturbance 
formulation. This allows the geometry of the configuration to be greatly 
simplified and the surface boundary condition to be applied on a mean approxi- 
mate surface. Accurate prediction of such flow fields, however, requires the 
use of the full-potential foraulation and the generation of a suitable 
surface-fitted grid. Because each aircraft component (wing, nacelle, fuse- 
lage) requires, in general, a grid system that is usually incompatible with 
the grid systems of the other components, the generation of a single surface- 
fitted grid for the entire configuration is a difficult task.. In such a 
global grid, control of grid point distribution, skewness, and clustering 
would be difficult to achieve. 


In the present investigation, an alternate approach is employed in com- 
puting the transonic flow field about wing/nacelle configurations. The 
approach taken in developing this multicomponent algorithm uses a component- 
adaptive grid embedding scheme in which the global computational grid is 
composed of a series of overlapped component grids, where each component grid 
is optimized for 3 particular geometry such as the wing or nacelle. The AF2 
approximate factorization algorithm 1- is used to determine the conservative 
full-potential equation solution on each component grid with trivariate inter- 
polation being used to transfer property information between the component 
grids. The AF2 algorithm has been applied to the computation of two- 
dimensional transonic airfoil flows by Holst^ and to the computation of three- 
dimensional transonic wing flows by Holst and Thomas -1 . Increases in conver- 
gence speed by factors of 4 to 7 have been realized using the AF2 scheme 
Instead of using the standard transonic relaxation scheme, successive-line- 
over-relaxation. 


Using the component-adaptive grid approach requires that the grid 
generation and flow solution algorithms be available for each component in the 
configuration. For determining the solution for wing/ nacelle . onf igurations, 
this then requires that an isolated wing algorithm and an isolated nacelle 
algorithm be available to be linked to form the combined multicomponent algo- 
rithm. In the present investigation, the NASA-Ames GRGEN3 wing grid gene- 
ration Fortran program and the TWING transonic wing flow analysis Fortran 
program -1 are employed in computing the wing component flow field. To deter- 
mine the nacelle component flow field, the NGRIDA nacelle grid generation 
program and the NACELLE nacelle flow analysis program are employed. The 
NGRIDA and NACELLE Fortran programs were developed as part of the present 
study and their use for the computation of three-dimensional transonic flows 
about isolated nacelle/ inlet configurations is documented in Reference 4. The 
TWING and NACELLE programs have been combined to produce a wing/nacelle multi- 
component flow analysis algorithm called TVJN. An additional Fortran program, 
called LK, has been developed which interfaces the component grids. 

Preliminary results employing this analysis method are reported in Reference 
5. 
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SECTION II 


COMPUTATIONAL METHOD 


1. INTRODUCTION 

Both the wing and nacelle body-fitted conputatlonal grids are determined 
using numerical mesh generation techniques. The GRGEN3 and NGRIDA Fortran 
computer programs determine the wing and nacelle component grids, 
respectively. After the component grids arc established, grid interfacing is 
performed using the LK Fortran program. The global computational grid is then 
used by the TWN Fortran flow analysis program in determining the multi- 
component flow field. The efficient AF2 approximate factorization scheme is 
used by the TWN algorithm in solving the system of finite-difference 
equations . 

In this section, the computational methods used to generate the 
curvilinear component meshes and to obtain the inviscid flow solution are 
presented. Computed results and comparison with experimental data are 
presented to illustrate application of the analysis. 

2. COMPONENT GRID TOPOLOGY AND GENERATION 

The global computational grid is composed of the wing and nacelle 
component grids, each of which is generated separately. Each component grid 
is body-fitted and is determined using a numerical grid-generation algorithm. 

2.1 WING COMPONENT GRID GENERATION 

The wing flow field solution is determined on a three-dimensional 
body-fitted curvilinear mesh. The wing computational mesh is obtained using 
two-dimensional numerical grid generation techniques for a series of spanwise 
stations which are distributed along the wing span starting at the root 
symmetry plane and progressing past the wing tip into the far field. 

Figure la illustrates the O-type body-fitted grid for an isolated wing. 
The base Cartesian coordinates are denoted by x, y, and z. The conputatlonal 
curvilinear coordinates are denoted by £, n , and C,. The £ coordinate is in 
the wraparound direction, initiates at" the wing upper surface trailing edge, 
and terminates at the wing lower surface trailing edge. The n coordinate is 
in the spanwise direction, initiates at the root symmetry plane, and 
terminates at the free-stream sidewall boundary. The S coordinate is in the 
radial or normal direction, initiates at the outer computational boundary, and 
terminates at the wing surface. 

The wing surface grid points are clustered near the wing leading and 
trailing edges using either a conformal mapping distribution or a geometric 
stretching function expressed in terms of the arc length measured from the 
wing edge points. The outer boundary points are distributed at equal angular 
increments around the body. 

Once the wing surface and outer boundary me-h point distributions have 
been determined for a given spanwise station, the interior mesh point 
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locations are deterained by iteratively solving the two coupled Laplace 
equations 


5 + £ 

XX zz 


0 


( 1 ) 


e + c 

XX *zz 


( 2 ) 


These equations are transformed to and solved in the computational domain 
where £ and c are the independent variables and x and z are the dependent 
variables. An alternating-direction-implicit (ADI) relaxation algorithm is 
used to solve the governing finite-difference equations^. This then 
establishes the grid point coordinates at each spanwise station used in 
defining the wing geometry. Clustering of the spanwise stations is performed 
in the vicinity of the wing root, tip, and span location of the nacelle. 

The GRGEN3 Fortran computer program^ is used to generate the wing 
component grid. For the case of a wing with no section variation, only the 
section geometry for the root defining station needs to be entered. In this 
case, the grid coordinates are found using the numerical relaxation technique 
for the root and tip stations. The grid coordinates for the remaining wing 
spanwise stations are found by linearly interpolating between the root and tip 
stations. Off the wing tip, the grid is generated for wing sections of zero 
thickness. For arbitrary wing configurations, taper, twist, thickness, and 
sweep variations can be accounted for by specifying the wing section geometry 
at a number of spanwise stations. 

Figure lb illustrates the rectangular computational domain which 
corresponds to the physical domain of Figure la. The physical space 
boundaries transform into boundaries in computational space. This allows for 
accurate and straightforward boundary condition implementation. The mesh 
points are equally spaced in computational space thereby permitting the use of 
standard finite-difference formulae in the flow solution analysis. 

To use the isolated wing grid in the present multicomponent flow field 
analysis, selected interior mesh points are deleted from the grid in order to 
accommodate the nacelle geometry. Figure 2 illustrates a view of the wing 
grid at a span station which contains the nacelle. The aperture in the 
interior of the grid in physical space creates an irregularly shaped void in 
the wing computational space domain. 

2.2 NACELLE COMPONENT GRID GENERATION 


The nacelle flow field solution is determined on a three-dimensional 
body-fitted curvilinear computational mesh. The nacelle computational mesh is 
obtained by using two-dimensional numerical grid generation techniques for a 
series of meridional planes which are splayed circumferentially around the 
body (a meridional plane is a plane containing the longitudinal axis of the 
inlet) . 

Figure 3 illustrates the C-type nacelle/ inlet body-fitted grid for 3 
nacelle extended in the downstream direction. The base Cartesian coordinates 
are denoted by x, y, and z. Figure 3a shows the meridional plane grid 
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Figure 2. Side view of wing component grid at a span 
station containing the nacelle. 
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topology, whereas Figure 3b Illustrates the grid as viewed along the 
longitudinal axis of the inlet which ia the local x axis. The conputational 
curvilinear coordinates are denoted by £ , n , and C . The £ coordinate is in 
the wraparound direction, initiates at the external outflow surface, and 
terainatea at the compressor face outflow surface. The 0 coordinate is in the 
circumferential direction, initiates at the 0*0 meridional plane, and 
terainatea at tha 0 *2ir neridional plane. The C coordinate is in the radial 
direction, initiates at the outer computational boundary (or centerline 
boundary), and terminates at the body surface. The 0*0 and 0*2^ 
meridional planes are coincident. If a pylon were present, then the pylon 
geometry would be contained between these two meridional stations. 

The surface grid points are clustered in the region of the nacelle hilite 
(leading edge point of the nacelle). The clustering is achieved by using a 
geometric stretching function which is expressed in terms of the arc length 
measured along the body from the hilite point. The outer computational 
boundary point distribution is determined by using either an angular 
distribution function expressed in terms of an angle measured about the 
nacelle hilite, or by using an arc length distribution along the outer 
boundary. 

Once the surface and outer boundary mesh point distributions have been 
determined, the interior mesh point locations are computed using the NASA-Anes 
GRAPE** algorithm. The GRAPE algorithm determines the interior field point 
coordinates by iteratively solving two coupled Poisson equations. For 
axisymaetric geometries, this two-dimensional grid generation procedure is 
applied for only one meridional plane, and the grid point locations on the 
remaining meridional planes are found by using simple reflection techniques. 
For asymmetric geometries, the two-dimensional grid generation algorithm can 
be applied for each neridional plane. 

The GRAPE two-dimensional grid generation algorithm is contained within 
the NGRIDA nacelle grid generation program and is used solely to obtain the 
interior grid point locations. The boundary point locations are computed in 
NCR IDA and are supplied to the GRAPE subalgorithm to be used as boundary 
conditions. After the interior field point coordinates have been determined, 
point reordering, and grid translation, scaling, and reflection are performed 
in NGRIDA. The detailed theory of the NGRIDA nacelle grid generation 
algorithm is discussed in Reference 4. 

Figure 4 illustrates the rectangular conputational domain which 
corresponds to the physical domain shown in Figure 3. The physical space 
boundaries transform into boundaries in computational space. 

Figure 5 illustrates a typical meridional plane interface grid. The 
inlet centerline, which is the x axis, represents a singularity in the three- 
dimensional grid mapping. The computational grid boundary is displaced a 
small distance away from the x axis. An extrapolation and averaging 
technique, described later, is used to obtain flow properties on the 
centerline. 

3. GOVERNING EQUATIONS FOR THE INVISCID FLOW 

The inviscid flow gas dynamic model is based on the assumption of steady 
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Figure 4. 


Nacelle computational domain. 
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potential flow which requires chat the flow be both irrotatlonal and 
l3cntroplc. The governing equations for steady three-dimensional potential 
flow are given by • \ / \ / > 



(3) 


(A) 


where 5 , H , and C denote the system of curvilinear coordinates, U, V, and W 
are the contravariant velocity components in the 5 , n , and Q directions,, 
respectively, p is the density, J is the Jacobian of transformation from the 
Cartesian coordinate system (x, y, z) to the general curvilinear coordinate 
system (S, n, C), 0 is the velocity potential function, and Y is the specific 
heat ratio. The density and contravariant velocity components are normalized 
by the stagnation density and critical sonic speed, respectively. 


Equation (3) is the full-potential equation in strong conservation law 
form. It expresses mass continuity for steady three-dimensional flows. 
Equation (A) expresses entropy conservation and is used to compute the density 
given the velocity potential field. 


The contravariant velocity components can be expressed in terms of the 
Cartesian velocity components as 
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where u, v, and w denote the nondimensional velocity components along the x, 
y, and z Cartesian coordinate axes, respectively. The contravariant velocity 
components can be also expressed in terms of the velocity potential 2 as 
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V 5 + Vn + V< 


( 8 ) 


v « + Vn + V; 


(9) 


M “ V? + Vn + V< 


( 10 ) 


The metric parameters (i * 1 to 6) and the Jacobian of transformation 
J are given by 
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The following notric rolatlono are also required In defining the napping from 
physical space to coaputatlonal space: 
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The metric parameters are obtained numerically using standard second- 
order or fourth-order accurate finite-difference formulae to compute 
derivatives of the form x r , x , x, , y r , etc. Then using the metric relations 
given by equations (13), ?he Inverse quantities 5 , n etc., are determined. 
Substitution of these values into equations (LI) and (12) yields the A^ and 
Jacobian J. 

4. NUMERICAL ALGORITHM 

4.1 BASIC FINITE-DIFFERENCE SCHEME 

The present numerical algorithm is based on the finite-difference 
formulation used by Holst and Thomas^ in computing transonic wing flows. In 
this algorithm, the full-potential equation residual is given by 

L$ . - 5 (rj) + <5 (i“ ) + 

i,j,k C i+h,j,k n J i.j-Hj.k 




( 14 ) 


i.j ,k+*j 


where Lb, , ^ denotes the residual operator, and i, j, and k denote the grid 
point indices in the E, (wraparound), n (circumferential or spanwise), and 5 
(radial) directions, respectively. The magnitude of the residual LO^ . ^ 
approaches zero as convergence is attained. The operators 6 r( ) , ( J ), and 0 , 

( ) are first-order accurate backward finite-difference operators (applied at ^ 
midpoints) for the £, h, and C directions, respectively. The terms p, p, and p 
are upwind-biased density coefficients given by expressions of the form 


p i+4.J,k ■ ta-v)pl i+4,j,k 


+ V i+ 1 3,j.k P i+4+r,J,k 


( 15 ) 


where r denotes an upwind point along the q direction, and v is an artificial 
viscosity coefficient. In equation (15), the physical density 0 is computed 
using equation (4) with central differences used for determining the 
derivatives of $ . The artificial viscosity coefficient v is given by 


> 

1 


12 
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v -° I? “ij.k* 1 


V -° ^t.j.1 


" l > IF 1 


( 16 ) 


where M is Che local Mach number, and C is a user-specified constant. The 
artificial viscosity coefficient C typically ranges from 1.0 to 2.0, with the 
larger values producing greater upuinding. Expressions similar to equation 
(15) hold for 7 and p which effect upwinding in the n and Q directions, 
respectively. 

The finite-difference equations are solved using the AF2 approximate 
factorization scheme which has proved to be significantly more efficient than 
successive-line-overrelaxation schemes^. The AF2 algorithm is written in a 
three-step form as: 


First-Step: 

(a - “ ~6 A 6 ) g? , ■* ctu)U" . . 

\ n j n i,J i.J.k 

+ \+l f i,j,k+l 


Second-Step : 

(A, I fl. 6. - - 5.A.5.) f? . . » g" 

* + <; s a ^ ; i,J ,* 5 i,j 


(17) 


(18) 


Third-Step: 


(a + <5 C )C 


n 

i.j.k 


.n 

t , . , 

i.J 


(19) 


In equations (17) to (19), a is a time-like factorization parameter chosen to 
maintain stability and attain fast convergence, 3r is a factor that controls 
the amount of dissipation required in regions of supersonic flow, u is a 
relaxation factor, n is the iteration number, Ld>. . k is the mass residual 
[defined by equation (14)], f and g are intermediate functions which are 
obtained during the solution process, and Cj . k is C he potential function 
correction given by 
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are defined by 



A l " <pA l /J) i-‘i,j,k 


(21) 

A j * (pA 2 /J) i,J-4,k 


(22) 

\ " (pA 3 /J) i,j,k-4 


(23) 


In Steps 1 and 2 the g and f functions are obtained by solving a 
tridiagonal systen of equations while in Stef, ; the correction ^ Is 

obtained by solving a bidiagonal systen of equations. 

The iterative relaxation procedure can be viewed as a tine-marching 
integration algorithm in pseudotime. Using this analogy, the factori- 
zation parameter a appearing in equations (17), (13), and (19) can be regarded 
as the reciprocal of the marching pseudotine step. To eliminate all 
components of the error frequency spectrum, it is generally preferable to 
employ a variable time step which sequentially varies with iteration number. 
The small values of a are particulary effective in reducing the low frequency 
errors, while the large values of ct are effective in reducing the high 
frequency errors. 


The factorization parameter a is computed using a repeating '-^uence 
which is based on iteration number'-. The ct sequence that is empj.oyed is given 
by 


k-1 

“k " °h (a l /a h } M_1 


(k - 1,2, — ,M) 


(24) 


where a ^ is the lower limit of a, is the upper limit, and is the value 

of a for the kth element of the sequence. For all cases presented herein, the 
number of elements in the sequence M was equal to 3. The optimum values of 
and are generally determined by numerical experiment. 

4.2 BOUNDARY CONDITIONS 

At the surface of the wing or nacelle, the contravariant velocity W in 
the 4-coordinate direction is equated to zero as this satisfies the flow 
tangency condition. At a body surface, the derivative tern $ is calculated 
from • ^ 


= - (A.5. + ) /A 

„ 5 c 6 n 3 


(25) 
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where the derivativeo and ^ are found using second-order differencing in 
the C"conatfint boundary surface. 

At tha conpreosor face outflow surface inside the inlet, the 
contravariont velocity U Q in the ^-coordinate direction is specified. Given 
the local Cartesian velocity components u, v, and w, and the metric quantities 
£ x , £ y , and S z , U Q can be determined using equation (5). Ac the compressor 
face, the flow is assumed to be uniform and in the axial (x) direction. The 
local flow velocity and density are fixed by specification of the required 
engine mass flow race. Once U Q has been determined, the derivative term is 
calculated from 

*5 = < U o - Vn ■ V; )/A i (26 > 


where 9 and 5 . are found using second-order differencing in the ^-constant 
boundary surface. 

For the nacelle component algorithm, the grid centerline boundary is 
offset a small distance from the centerline to avoid the mapping singularity 
at the centerline. To determine property values on the centerline, an extra- 
polation and averaging procedure is employee, which is fully discussed in 
Reference 4. 

At the wing symmetry plane boundary, the contravariant velocity V in the 
q -coordinate direction is equated to zero. At Che wing free-streaa sidewall 
boundary, the initial free-stream distribution of 0 is held fixed as is the <? 
distribution for the wing grid outer boundary for non-lifting cases. For 
lifting calculations, the 9 distribution on the wing grid outer boundary is 
updated using a compressible vortex solution as this improves convergence 
speed 3 . The wing circulation is computed at the end of each wing algorithm 
iteration by calculating the potential jump along the vortex sheet. 

■ The nacelle component grid outer boundary constitutes the grid overlap 
region outer boundary. The grid overlap region inner boundary is constituted 
by the irregularly shaped boundary which is obtained from deleting selected 
interior field points from the wing component grid in order to accommodate the 
nacelle geometry (see Figure 2). For the overlap region inner and outer 
boundaries, the velocity potential 9 or the velocity potential normal 
derivative <> n can be specified. Both options were tried in two dimensions 2 
and no significant differences in the final results were observed. However, 
specifying the velocity potential which is a scalar quantity is easier to 
implement and requires less computational effort as compared to specifying the 
velocity potential normal derivative. 

The portion of the overlap region inner boundary adjacent to the engine 
exhaust plane at the aft end of the nacelle requires special treatment. A 
uniform axial velocity is specified at the nacelle aft station and 3 
three-point extrapolated difference is used to provide the required values of 
the potential function at the wing grid points immediately downstream of this 
station. 
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In addition to (specifying $ on the overlap region boundaries, it io 
nececoary to opocify boundary conditiono for the f and g interoediate 
corrections for uea in the wing component algorithm. For the f Intermediate 
correction in the 5-wraparound coordinate direction, the Dirichlet boundary 
condition f ■ 0 is imposed on the overlap region inner boundary in performing 
the wing component flow field calculations. For the g intermediate correction 
in the n-opanvioo coordinate direction, the Neumann boundary condition g ■ 0 
is imposed on the overlap region inner boundary. Both the f and g corrections 
approach zero as the iterative procedure progresses and convergence is 
attained. 

For the nacelle component algorithm, the Dirichlet boundary condition f ■ 
0 is imposed at the external outflow surface of the nacelle grid in solving 
the system of factorized difference equations in the q-wraparound coordinate 
direction. Periodic differencing is employed in the n-circumferential 
direction. 


4.3 ITERATION PROCEDURE 

The computation is started by initializing the entire velocity potential 
field. For the external flow around the wing and nacelle cowling, the 
initialization of the potential field is performed using free-stream velocity 
components. For the internal flow initialization inside the inlet, the local 
velocity is assumed to be axial and is computed using a Mach number which is 
determined from the implicit relation 
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where M and M C £ denote the local and compressor face Mach numbers, 
respectively, and A and A c denote the local and compressor face flow areas, 
respectively. Equation (27) was obtained using one-dimensional gas dynamic 
formulae and simply expresses mass continuity for the captured stream tube. 

To determine the internal potential field, a trapezoidal rule integration is 
used. The compressor face Mach number M c £ is fixed by specification of the 
required engine mass flow rate. 

Once the entire potential field has been initialized, the ving algorithm 
is then executed for a specified number of iterations (typically 4? -50) 
holding $ constant on the overlap inner boundary. The nacelle algorithm is 
then executed for a specified number of iterations (typically 40' 50) using 
overlap outer boundary conditions as determined from the wing solution. \<. 
this stage, the wing algorithm is again executed using updated overlap inrv.r 
boundary conditions obtained from the nacelle solution. This process is 
repeated until overall convergence is achieved. Typically 4 to 5 cycles are 
required for convergence with 40 to 50 iterations being performed in both the 
wing and nacelle component algorithms per cycle. 

4.4 INTERGRID PROPERTY TRANSFER 

Multivariate interpolation is employed to transfer property infor-cation 
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between tha component grids. To effect intergrid property transfer for the 
wing/nacallo configuration requires employing a trivariate interpolation 
polynomial oincc tho component grids do not have a common coordinate surface. 
To reduce tha coding logic, a trivariate polynomial based on a linear Taylor 
series expansion was employed. Interpolation is performed in computational 
(£,0»C) space using a polynomial of the form 


f(5,n,0 - «„ + <5 -5 0 > + f <n-n 0 > + f. «-?„> <«> 

o 

where f(5,n,C ) denotes the interpolated flow property at the point ( C, 0 , s ), 
and the subscript o denotes a selected interpolation base point in 
computational space. The derivatives appearing in equation (28) are evaluated 

usine standard second-order accurate difference formulae. Point ( S „ Z ) 

° o * O 9 \ 

is selected as the mesh point closest to the desired interpolation point (^, h, ->) • 


The use of the trivariate interpolation polynomial given by equation (28) 
in some cases produced oscillations in the interpolated values of the velocity 
potential. The nacelle component grid lies in a region where the wing 
component grid is relatively coarse, thereby creating a disparity in the 
average cell size between the two grids in the overlap region where the 
interpolation is performed. It was found that the best results were obtained 
by employing a relatively fine grid for the wing and a relatively coarse grid 
for the nacelle. To ensure a monoconic behavior for the interpolated overlap 
boundary potential values, explicit smoothing is performed in the streamwise 
direction after each cycle of the wing/nacelle solution. 

5. NUMERICAL RESULTS 

Selected numerical results are now presented to Illustrate application of 
the analysis. 

The first set of computed results are lor a hypothetical wing/nacelle 
configuration. The wing geometry is represented by a NACA 0012 wing which has 
an aspect ratio of 6.0, zero twist, and which is unswept and untapered. The 
nacelle geometry is represented by the Lockheed-Georgia GELAC1 axisymmetric 
nacelle/ inlet configuration. The nacelle is located at the 50% span station 
and has its axis parallel to the wing chord. 

The wing geometry is the default geometry contained within the GRGEN3 
wing grid generation code. The nacelle geometry is the default geometry 
contained within the N'GRIDA nacelle grid generation code. The results of 
isolated nacelle Clow field computations for the GELAC1 nacelle are presented 
in Reference 4. 

The grid topologies used for this wing/nacelle multicomponent flow field 
computation are illustrated In Figures 6, 7, and 8. Figure 6 shows a frontal 
view of the composite overlapped grid system. The outer wing grid is 
generated for a series of spanwise stations. The span stations are clustered 
near the location of the nacelle. The inner nacelle grid is generated for a 
series of meridional planes spaced in equal angular increments around the 
body. A side view of the composite grid system at the span station coincident 
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Figure 7. Side view of NACA0012-wing/GELACl-nacelie 
configuration global grid. 
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with tha oscallc oycnatry plane is given in Figure 7. The overlap region 
between the wing O-grid and tha nacelle C-grid is clearly illustrated in this 
figure. Tha nacelle interface grid is illustrated in isolation in Figure 8. 
Tho grid generation procedure used to obtain a nacelle interface grid i3 
alnost identical to that used to obtain an isolated nacelle grid except Chat 
for an interface grid the outer computational boundary is placed close to the 
body and typically the number of radial stations is reduced. The wing 
component grid for this case consisted of 101 wraparound stations, 22 spanwise 
stations, and 30 nornal stations. The nacelle coaponent grid consisted of 67 
wraparound stations, 13 circuaferential stations, and 5 radial stations. 

Figure 9 illustrates the conputed surface solution for this wing/nacelle 
configuration at a free-strean Mach nuaber M^of 0.63 and at an angle of 
attack a of 2.0 degrees. Shown in this figure is the surface pressure 
coefficient c versus the fractional distance along the chord x/c (distance/ 
chord length) p for the 50% span station (H » 0.5). Results are presented for 
the wing/nacelle multicomponent flow calculation and for an isolated wing 
calculation. The wing/nacelle calculation was performed assuming that the 
inlet capture ratio was 0.7 and that the ratio of the engine exhaust plane 
velocity to the free-stream velocity was 2.25. The effect of including the 
nacelle is to reduce the suction level on the upper wing surface and to 
increase it on the lower surface near the aft portion of the nacelle (wing 
leading edge). The high engine exhaust Jet velocity increases the flow 
velocity near the lower wing surface, thereby reducing the circulation and 
hence the lift. The midspan wing section lift coefficient for the isolated 
wing calculation was 0.287. For the wing/nacelle calculation, the midspan 
section lift coefficient was reduced to 0.24. These results agree 
qualitatively with what is observed experimentally®. The wing/nacelle mutual 
interference calculation required four cycles through both the wing and 
nacelle component algorithms to achieve convergence which corresponds to a 
three order of magnitude reduction in the maximum residual. The computation 
required approximately 5 hours of CPU execution time on the DEC VAX-11/780 
computer. Approximately one million words of core storage were required for 
the computation. 

Figure 10 illustrates the computed midspan wing pressure distribution for 
the XACA0012-wing/GELACl-inlet configuration for the conditions of a 
free-stream Mach number of 0.75 and 1 degree angle of attack. For this 
computation, the inlet capture ratio was maintained at 0.7 but the ratio of 
the engine exit plane velocity to the free stream velocity was now equated to 
1.15. Again, results are presented for an isolated wing solution and a 
w.tng/nacelle mutual interference solution. In this case, the wing upper 
surface pressure distribution is slightly altered by the presence of the 
nacelle, however, the lower surface pressure distribution is markedly 
different. In this case, the lower surface flow undergoes a rapid expansion 
followed by a strong shock at the 40% chord location. The net effect is a 
reduction in the sectional lift. These results indicate that the sense and 
magnitude of the nacelle interference effects are strongly dependent on the 
exhaust jet velocity, free-stream conditions, and nacelle location. 

The next set of computed results are for a recently-designed 
Lockheed-Georgia wing/nacelle configuration. The wing geometry employs 
supercritical airfoil sections and is both swept and tapered. The nacelle 
geometry is again the GELAC1 axisymmetric nacelle/ inlet configuration. The 
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Figure 9. Wing pressure distribution for 

NACA0012-wing/GELACl-nacelle conf iguration 
at M *■ .63 and a = 2°. 


Figure 10. Wing pressure distribution for 

NACA001 2 -wing/OELACl -nacelle conf iguration 
at M * .75 and a *• 1° . 
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nacelle plane of aymmatry is again located at the wing aidspan ( n* 0.5). 
Flgureo li, 12, and 13 illustrate tho global coaputational grid for this wing/ 
nacelle configuration. Figure 11 chows a frontal view of the overlapped 
ccaponent grid system, whereas Figure 12 9hows a aide view at the span 
location of tho nacelle symaatry plane. Figure 13 illustrates a planfona view 
of the wing and nacelle. The wing ccaponent grid employed 101 wraparound 
stations, 22 spanwise stations, and 30 normal stations. The nacelle component 
grid employed 67 wraparound stations, 13 circumferential stations, and 5 
radial stations. 

Figures 14, 15, and 16 present the computed wing pressure distributions 
for this wing/nacelle configuration at the 463, 503, and 603 span stations, 
respectively. The free-stream Mach number and incidence used in the 
computations correspond to 0.77 and 2 degrees, respectively. On each of the 
three figures, the computed solution is shown for an isolated wing calculation 
and for wing/ nacelle calculations using exhaust jet to free stream velocity 
ratios (v^/v^) of 1.2 and 1.9. The inlet capture ratio was 0.7. The nacelle 
interference effect for each of the three spanwise locations was to reduce the 
suction level on the wing upper surface with this effect being more pronounced 
as the exhaust jet velocity was increased. As the upper surface suction level 
was reduced, the shock was moved slightly forward for each of the three span 
locations. The interference effect of the nacelle on the wing lower surface 
was to increase the suction levels near midchord. Again this effect was more 
pronounced at the higher exhaust jet velocity. On the lower surface near the 
wing leading edge, the surface pressure was higher for the wing/nacelle 
configuration than the isolated wing except at the inboard station where the 
pressures were comparable. For all three span stations illustrated, the 
interference effect of the nacelle decreased the sectional lift coefficients 
as compared to the isolated wing results. Moreover, the reduction in 
sectional lift became more pronounced as the exhaust jet velocity was 
increased. Overall, the computed effects of the nacelle interference with Che 
wing flow field agree qualitatively with experimental observations. 

The final set of results are for the Lockheed C-5A military transport 
aircraft^. Figures 17, 18, and 19 present a frontal view, side view, and 
planform view of the overlapped grid system used for this computation. Figure 
20 illustrates the C-5A nacelle interface grid. The contours of the nacelle 
used in the computation represent a flow-through nacelle which is used to 
simulate the actual turbofan flight ometry which is illustrated in Figure 
21. The wing component grid employee 101 wraparound stations, 22 spanwise 
stations, and 30 normal stations. The nacelle component grid employed 67 
wraparound stations, 13 circumferential stations, and 5 radial stations. 

Figures 22 and 23 illustrate the computed and experimental wing pressure 
distributions at the 393 and 463 span stations, respectively. These results 
are for a free-stream Mach number of 0.775 and zero incidence. The nacelle 
capture ratio was assumed to be 0.55 and the engine exit plane exhaust 
velocity was taken as 1.125 times that of the free-stream velocity, figure 22 
presents computed and experimental results for both an isolated wing and a 
wing/nacelle configuration. Figure 23 presents results for the wing/nacelle 
case only. Referring to the experimental results shown on Figure 22, the 
effect of adding the pylon and nacelle is to decrease the peak suction level 
on the wing upper surface. This trend is predicted by the analysis 
qualitatively. The analysis, however, overpredicts the suction level on the 
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(Vj/V^ = 1.2) 



Figure 14. Wing pressure distribution for Lockheed- 
Georgia wing/nacelle configuration at 
46% span. 
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Figure 15. Wing pressure distribution for Lockheed 
Georgia wing/nacelle configuration at 
50% span. 
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Figure 16. Wing pressure distribution for Lockheed- 
Georgia wing/nacelle configuration at 
60% span. 
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Figure 17. 


Front view of Lockheed C-5A wing/nacelle 
configuration global grid. 
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Figure 18. Side view of Lockheed C-5A wing/nacelie 
configuration global grid. 



32 







I 

I 


CRiGtfJAL FAC' F5 
OF POOR QUALlTV 



"3 

•U 

u 

CO 

C 

3 

Z) 


0 
u 
C3 

c 

< 

in 

1 

o 

•o 

o • 
CJ c 
pG o 

US -H 
U *J 

O 3 
►J w 
3 

4* G. 
O £ 
O 

3 O 

<u 

c 

> 

O TJ 
*U CJ 

•u co 

tO 3 


o 

CN 

<y 

t- 

3 

CO 

tu 


3 ^ 


ORIGINAL PAGE E3 
OF POOR QUALITY 



a 

u 

n 

c 

c 

a 

w* 

o 

•O 

J-i 

3 


I 

O 

•a 

V 

Zi 

u 

o 

►J 


V 
u 
rs 

'** 3 
O "3 
U 

3 a 
a .c 


> « 
3: 
a eo 
-a 


co 


a 

u 

3 

CO 



35 












OF POC'A 1 

'ppar ourfaca from approximately the 302 to 702 chord stations* This 
overprcdictlon io prceont both In the Isolated wing and wing/nacelle 
computations. Tha agreement between analyais and experiment for tho wing 
lower nurfaca to substantially better. The analyais overpredict3 the lower 
surface static pressure near the leading edge. 

The difference between the results of the analysis and experiment on tho 
wing upper surface are attributed mainly to using conservative differencing 
without employing a boundary layer correction and to wind tunnel Interference 
effects. The usa of conservative differencing generally predicts stronger 
shocks than nonconservative differencing. However, invoking a boundary layer 
correction will generally correct the computed shock strength. The angle of 
attack used in the analysis was the same as that of the experiment. The 
experimental wing lift coefficient was not available^. Consequently, it was 
not possible to determine the angle of attack to he used in the analysis that 
would match the measured lift coefficient and then to compare the pressure 
distributions at the corrected incidence'. The difference between the results 
of the analysis and experiment on the wing lower surface near the leading edge 
are attributed mainly to the absence of pylon modeling in the analysis. 

Figure 24 compares the results of tho analysis and experiment for the 
C-5A configuration at the 602 span station. This case corresponds to a 
froe-stream Mach number of 0.7, a 2 degree angle of attack, an inlet capture 
ratio at 0.55, and exhaust jet to free-stream velocity ratio of 1.2. Results 
are presented tor an isolated wing and a utng/nacelle configuration. In this 
case, the agreement between analysis and experiment is encouraging on the wing 
upper surface for the wing/nacalle configuration. The -agreement on the wing 
lower surface is also good in the region from the wing aidchord aft. 

An effort was made to include swept pylon effects into the analysis. A 
pylon of infinitesimal thickness was positioned at the spanwise station 
containing the nacelle symmetry plane. A flow tangency condition was imposed 
on the pylon surface enforcing the contravariant velocity V in the n-direction 
to be taro on both sides of the pylon. This condition was imposed in both the 
wing and nacelle component algorithms. Difficulties were encountered, 
however, in achieving converged solutions with the pylon model included. This 
was felt to be due to the fact that the srid was not aligned with the swept 
pylon leading edge and to the fact that with wing sweep a spanwise velocity 
exists which causes the pylon to be loaded thereby requiring the inclusion of 
a vortex sheet which emanates from the pylon trailing edge. Time constraints 
in the present investigation did not allow these difficulties to be completely 
resolved. 
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Figure 24. Wing pressure distribution for Lockheed C-5A 
wing/nacelle configuration at 60" span. 
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SUBROUTINE DESCRIPTIONS 


l. INTRODUCTION 


In this section, a brief description is given of the function of each 
subroutine in GRGEN3 and NGRIDA grid generation Fortran programs, the LK grid 
interface Fortran prograta, and the TUN flow solution Fortran progran. This 
information supplements the information available in the form of comment 
statements within the respective programs. 

2. SUBROUTINE DESCRIPTIONS FOR THE GRGEN3 WING GRID GENERATION PROGRAM 

GRGEN3 . This program routine is the main control routine in the GRGEN3 wing 
grid generation program. The program execution is initiated by 
entering NAMELIST GRID3D which defines the parameters needed to 
determine the three-dimensional body-fitted wing grid. Routine 
GRGEN3 then calls SUBROUTINE YDIST to establish the grid plane dis- 
tribution in the spanwise direction. SUBROUTINE GRGEN is then called 
to determine the grid point distribution for each spanwise station 
that is used to define the wing geometry. The finished grid coordi- 
nate data are then loaded onto disk files TAPE 9 and TAPE 10 for 
input to the LK grid interface program. 

YDIST. This subroutine is used to establish the spanwise grid coordinate 
distribution. The spanwise grid planes are usually clustered near 
the wing root and wing tip and near the span location of the nacelle. 
SUBROUTINE YDIST is called from routine GRGEN3. 


GRGEN . This subroutine is used to determine the grid point distribution on a 
given spanwise station. NAMELIST GRIDIN is first read to enter the 
parameters required in the ADI relaxation scheme. SUBROUTINE GEMPAC 
is then called to determine the coordinates of the local wing 
section. At this stage, SUBROUTINE INNER is called to ascertain the 
surface grid point distribution by clustering points near the leading 
and trailing edges of the wing. SUBROUTINE OUTER is then called to 
determine the outer boundary grid point locations. Finally, SUB- 
ROUTINE GRGEN calls SUBROUTINE ADI to determine the interior field 
point locations by iteratively solving two coupled Laplace equations. 
SUBROUTINE GRGEN is called from routine GRGEN 3 . 

GEMPAC . This subroutine determines the wing section coordinates for each span 
station Chat is used to define the wing geometry. Seven basic geom- 
etry definitions are contained within SUBROUTINE GEMPAC and are con- 
trolled by the input parameter I0PT. SUBROUTINE GEMPAC is called 
from SUBROUTINE GRGEN. 

INNER ■ This subroutine is used to determine the grid point distribution on 

the wing surface for each spanwise station at which the wing geometry 
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io defined. The surface grid point clustering la controlled by the 
BI1JH input parameter. SUBROUTINE INNER la called from SUBROUTINE 
GRGEH. 

OUTER. Thio subroutine la used to determine the grid point distribution on 
the circular outer computational boundary for each spanwise station 
at which the wing geometry 13 defined. The circular outer boundary 
is centered at the wing section leading edge. SUBROUTINE OUTER Is 
called froa SUBROUTINE GRGEN. 

GOUT. This subroutine Is used to print the computed grid point di3tri- 

’ button. The amount of printout is controlled by the input parameters 
N0UT5, IINCR, and JINCR. 

TRIB . This subroutine is used to solve a tridiagonal system of simultaneous 
linear equations with fixed boundary conditions. SUBROUTINE TRIB is 
called from SUBROUTINE ADI. 

TRIP . This subroutine is used to solve a tridiagonal system of simultaneous 
linear equations with periodic boundary conditions. SUBROUTINE TRIP 
is called from SUBROUTINE ADI. 

■ADI . This subroutine is used to determine che final interior grid point 

distribution by solving two coupled Laplace equations using an alter- 
nating-direction-implicit iteration algorithm. SUBROUTINE ADI is 
called from SUBROUTINE GRGEN. 

CSPLIN. This subroutine is used to perform cubic spline interpolation to 

obtain the final distribution of surface grid points at each spanwise 
station where the wing geometry is defined. SUBROUTINE CSPLIN is 
called froa SUBROUTINE INNER. 

IN'TERP . This subroutine is used to perform linear interpolation for the grid 
point coordinates between the spanwise stations at which the wing 
geometry is defined. SUBROUTINE INTERP is called from routine 
GRGEN3. 

3. SUBROUTINE DESCRIPTIONS FOR THE NCR IDA NACELLE GRID GENERATION PROGRAM 

NGRIDA. This program routine is the main control routine in the NGRIDA grid 

generation program. Routine NGRIDA first calls SUBROUTINE NINPUT for 
data input, parameter initialization, and printing of preliminary 
information. SUBROUTINE GENER is then called to control the actual 
numerical generation of the surface-fitted computational grid. 

NINPUT . This subroutine is used to enter the input data, perform parameter 
initialization, and print preliminary information. The input data 
are entered by specification of the four NAMELISTS LIST1, LIST2, 
LIST3, and LIST4. All input parameters have specified default 
values, and a default axisyrametric nacelle/ inlet geometry d3ta set 
has been loaded onto DATA statements within the program. At the end 
of SUBROUTINE NINPUT, the input namelist data are printed to initiate 
the computer output. SUBROUTINE NINPUT is called from routine 
NGRIDA. 
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CUBICS . Thio subroutine determines the coefficients for a set of cubic spline 
polynomials uhich are used to interpolate for tho surface geometry 
given the axial positions and corresponding surface radii for N tabu- 
lar data points. The use of cubic splines ensures continuity of 
surface radius, slope, and curvature for the interpolated geometry. 
Tho oplino coefficients are found in terns of the surface curvature 
values at each of tho supplied tabular data points. A sparse systen 
of 3inultaneous linear equations, uhich deternine3 the curvature 
values, is solved by calling SUBROUTINE GELG. SUBROUTINE CUBICS is 
called from SUBROUTINE GENER. 


GENER. This subroutine contains the main control logic for nunerically 

generating the nacelle curvilinear computational grid. The initial 
step in GENER is to generate the cubic spline polynomial coefficients 
for both the external and Internal nacelle surfaces by calling SUB- 
ROUTINE CUBICS. The cubic spline polynomials are used to interpolate 
for the surface geometry. At this stage, the arc length measured 
from the nacelle hilite is determined at each of the user supplied 
tabular data points that define the geometry. From a set of user- 
supplied stretching factors uhich effect clustering around the 
hilite, the arc length distributions for the final surface grid point 
locations are determined. From the surface grid point arc length 
distributions, the axial and radial coordinates of the final surface 
grid points are found by interpolation using the cubic spline poly- 
nomials. At this stage, the outer computational boundary grid point 
distribution is determined by one of two options. The first option 
determines the outer boundary point locations by using an angular 
distribution where the boundary points are located on rays spaced in 
equal angular increments around the nacelle hilite. The second 
option determines the outer boundary point locations by using an arc 
length distribution with the spacing between points being controlled 
by user-supplied stretching parameters. At this stage, a file is 
written for the GRAPE algorithm which determines the interior field 
point locations by iteratively solving two coupled Poisson equations. 
After the GRAPE algorithm has been executed, conversion of the data 
format and plotting of the finished grid is performed by calling 
SUBROUTINE CONVERT. At this st3ge, a reflection of the grid point 
coordinates is performed to determine the point distributions on the 
remaining meridional planes. Finally, certain transition point in- 
dices are determined and the grid point distribution is printed and 
also stored on a disk file for input to the flow solution program. 
SUBROUTINE GENER is called from routine N’GRIDA. 


CONVERT . This subroutine is used to convert the output from the GRAPE algo- 
rithm into a format which is acceptable by the TWN flow solution 
algorithm. The GRAPE-generated grid point coordinates are entered 
through a binary read of TAPE 7. The corners on the leading edge of 
the grid are then rounded and point reordering is performed to make 
the grid generator output consistent with the input format required 
by the TWN program. Translation and scaling of the grid are then 
performed. If desired, plotting of the finished grid is performed. 
The plotting routines use the ISSCO-DISSPLA 10 Fortran library of plot 
functions. If this subroutine library is not available, then the 
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plotting ooctlon of the code will have to be modified to be conpat- 
lblo with thQ uoor'o particular lnatallatlon. SUBROUTINE CONVERT Is 
called froa SUBROUTINE GENER. 

CELG. Thlo oubroutlna Is u»ed for solving a system of simultaneous linear 
equationo. Tho system la oolved using Gaussian elimination with 
complete pivoting. Thlo subroutine is IBM Library SUBROUTINE GELG. 
SUBROUTINE GELG is called froa SUBROUTINE CUBICS. 

The remaining subroutines in the NGRIDA grid generation program com- 
prise the GRAPE algorithm. Cursory descriptions of the functions of these 
subroutines are presented herein. Detailed subroutine descriptions are given 
in Reference 6. 


GRAPE. 


CKSMTH . 


CSPLIN. 


IC. 


INCHK. 


INNER. 


INPUT . 


INTERP. 


OUTL.i. 


OUTPUT . 


RELAX. 


This subroutine is the main control routine in the GRAPE algorithm. 
SUBROUTINE GRAPE is called from SUBROUTINE GENER after the boundary 
grid point locations have been determined. SUBROUTINE GRAPE deter- 
mines the interior field point locations by calling in turn SUB- 
ROUTINES INPUT, INCHK, INNER, OUTER, SOLVE, and OUTPUT. The grid 
point locations are transmitted back to SUBROUTINE GENER through u3e 
of a scratch disk file. 

This GRAPE algorithm subroutine is used to check for smoothness in an 
input array by passing a parabola through the three nearest neighbors 
of a point and then observing the difference between the actual value 
and that interpolated from the parabola. 

This GRAPE algorithm subroutine is used to perform cubic spline curve 
fitting. 

This GRAPE algorithm subroutine is used to initialize the grid coord- 
inate and forcing function arrays prior to starting the iterative 
relaxation procedure to determine the interior grid point locations. 

This GRAPE algorithm subroutine is used to test for errors on select- 
ed input parameters for the GRAPE program. 

This GRAPE algorithm subroutine is used to locate grid points on the 
inner computational boundary. 

This GRAPE algorithm subroutine is used to enter the input data that 
is required for execution of the GRAPE program. 

This GPAPE algorithm subroutine is used to interpolate from the 
coarse grid solution to obtain initial conditions for the fine grid 
solution. 

This GRAPE algorithm subroutine is used to locate grid points on the 
outer computational boundary. 

This GRAPE algorithm subroutine is used to print the finished grid 
point coordinates. 

This GRAPE algorithm subroutine is used to solve the governing Pois- 
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son equations using a succesoive-line-over-relazation iteration 
ochcma < 

SOLVE. This GAPE algorithm subroutine is used to control the coarse and 
fins grid solutions of the governing Poisson equations. 

TRIB. This GEAFE algorithm subroutina is used to oolve a tridiagonal system 

““ of sinultaneous linear equations with fixed boundary conditions. 

TRIP. This CRAPE algorithm subroutine is U3ed to solve a tridiagonal systen 
of sinultaneous linear equations with periodic boundary conditions. 

4. SUBROUTINE DESCRIPTIONS FOR THE LK GRID INTERFACE PROGRAM 

LK. This program routine is the main control routine in the LK grid 

interface program. Routine LK enters the wing grid coordinate data 
by reading disk files TAPE 9 and TAPE 10 along with the nacelle grid 
coordinate data by reading disk file TAPE 14. From the parameters 
entered on NAMELIST NACEL, routine LK then scales the nacelle coord- 
inate data to be consistent with the wing reference length. The 
nacelle grid coordinates are then transformed into the wing coord- 
inate system. At this stage, the indices of the overlap region inner 
and outer boundaries are ascertained. SUBROUTINES FIX and MIX are 
then called to determine the interpolation grid point stencils for 
the points on the overlap boundaries. Finally, the interpolation 
point stencils and other variables required for interfacing the two 
grids are stored on the disk files TAPE 19 and TAPE 44. 

MIX. This subroutine is used to search the wing grid to ascertain those 

points that are nearest to a given point on the outer boundary of the 
nacelle grid. SUBROUTINE MIX also performs a search of the nacelle 
grid to ascertain those points that are nearest to a given point in 
the wing grid that i3 on the overlap region inner boundary. SUB- 
ROUTINE MIX is called from routine LK. 

FIX. This subroutine ascertains the interpolation point stencils used to 

determine the coefficients in a linear trivariate interpolation poly- 
nomial that is employed to transfer property information between the 
component grids. The interpolation polynomial s based on a first- 
order accurate trivariate Taylor series expansion. SUBROUTINE FIX is 
called from routine LK. 

5. SUBROUTINE DESCRIPTIONS FOR THE TWN FLOW SOLUTION PROGRAM 

TWN. This program routine is the main control routine In the TWN flow 

analysis program. Routine TWN calls SUBROUTINES INPUT and INIT to 
enter the input data and to perform the parameter initialization 
which are required for the computation of the wing component flow 
field. SUBROUTINES NINPUT and NINIT are then called to enter the 
input data and to perform the parameter initialization which are 
required for. the computation of the nacelle component flow field. 

The wing and nacelle flow field solutions are then computed in a 
cyclic iterative canner until overall convergence is achieved. SUB- 
ROUTINES INTW and INTO are called in the course of the computation to 
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INPUT. 


NINPUT. 


INIT . 


INTW. 


INTC.' 


NINIT. 
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interchange property data froa the wing and nacelle coaponent grids. 

This subroutine io used to enter the input data and perform initiali- 
zation of oolccted parameters which arc used in the coaputation of 
tha wing coaponent flow field. The input data arc entered by speci- 
fication of tha NAMELIST FLOUIN. The coordinate data for the wing 
coaponent grid are entered by reading disk files TAPE 9 and TAPE 10. 
SUBROUTINE INPUT is called from routine TUN. 

Thi3 subroutine is used to enter the input data and perfora initiali- 
zation of selected paraneter3 which are used in the computation of 
the nacelle component flow field. The input data are entered by 
specification of the NAMELIST FLOWC. SUBROUTINE NINPUT is called 
from routine TUN. 

This subroutine is used to initialize selected parameters employed in 
the computation of the wing component flow field. SUBROUTINE INIT 
reads disk file TAPE 19 to enter the grid overlap region coordinate 
data that are generated by the LK grid interface program. SUBROUTINE 
INIT then calls SUBROUTINES METRIX and WRITMS to compute and store, 
respectively, the metric parameters used in the wing grid mapping 
from physical space to computational space. Finally, the wing 
potential field is initialized using free-streaa velocity components. 
SUBROUTINE INIT is called from routine TUN. 

This subroutine is used to interpolate the velocity potential values 
from the nacelle flow field solution which are then used as boundary 
conditions on the overlap region inner boundary for an ensuing wing 
flow field calculation. SUBROUTINE INTW is called KCX tines after 
every ITMAX iterations performed in the nacelle grid. SUBROUTINE 
INTW is called from routine TUN. 

This subroutine is used to interpolate the velocity potential values 
from the wing flow field solution which are then used as boundary 
conditions on the overlap region outer boundary for an ensuing 
nacelle flow field calculation. SUBROUTINE INTC is called KCX times 
after every MAXIT Iterations performed in the wing grid. SUBROUTINE 
INTC is called froa routine TUN. 

This subroutine is used to initialize selected parameters employed in 
the computation of the nacelle component flow field. SUBROUTINE 
NINIT first computes functions of the specific heat ratio and bi- 
nomial expansion coefficients used in the computation of the physical 
density. SUBROUTINE METRIC is then called to determine the metric 
parameters used in the nacelle grid mapping from physical space to 
computational space. At this stage, the free-stream velocity magni- 
tude and Cartesian velocity components are computed along with the 
appropriate outflow boundary condition parameters. Finally, the 
potential field is initialized. The external flow field points are 
initialized using free-stream velocity components, whereas the inter- 
nal flow field points are initialized using velocities based on a 
local Mach number which itself is found from the area ratio and 
compressor face Mach number using one-dimensional gas dynamic formu- 
lae. SUBROUTINE NINIT is called from routine TUN. 
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Thin eubroutino la the main wing flow field integration control rou- 
tiaa end in used to apply tho AF2 approximate factorization ocheao to 
deteralne tho invincid wing flow solution. The AF2 schene is applied 
in an iterativo manner with the iteration sequence being terminated 
by either achieving convergence or by reaching a maximum perniaoible 
number of iterations. Per iteration, a three-step process is used to 
obtain the potential function correction. The first two steps calcu- 
late intermediate potential function corrections arising froa solving 
systems of factorized difference equations in the 5-wraparound and 
q-spanwise directions. The third step calculates the actual poten- 
tial correction by solving a system of factorized difference equ- 
ations for the 5 -radial coordinate. Steps l and 2 require sweeping 
the computational mesh from the body surface outwards. The third 
step sweeps from the outer computational boundary towards the body 
surface. SUBROUTINE SOLVE calls SUBROUTINES RO, ROCO, and RESID for 
calculation of the physical density coefficients, modified density 
coefficients, and mass residuals, respectively. SUBROUTINES TRIB and 
TRIP are called in the course of applying the AF2 scheme to solve 
systems of tridiagonal simultaneous linear equations with fixed and 
periodic boundary conditions, respectively. SUBROUTINE SOLVE is 
called by routine TWN. 

This subroutine is used to compute the metric parameters for all 
points in the wing component grid. The metric coefficients define 
the three-dimensional grid napping from physical space to computat- 
ional space. The metric coefficients are computed using fourth-order 
accurate finite-difference formulae. SUBROUTINE METRIX is called 
from SUBROUTINE INIT. 

This subroutine is used to compute and store the metric quantities 
that are used in computing the mass residuals at the wing surface. 
SUBROUTINE METCO is called from SUBROUTINE SOLVE. 

This subroutine is U3ed to calculate the physical density coeffici- 
ents at the 5-half points in the wing computational mesh. From the 
current potential function field, the derivatives of the velocity 
potential in the 5,0, and s directions are computed using second- 
order accurate finite-difference formulae. The potential function 
derivatives are then used in a form of the isentroplc energy equation 
to determine the density. At the body surface, the flow tangency 
condition is used in computing the density. At the external compu- 
tational boundary, the density is assumed to be equal to the free- 
stream density. SUBROUTINE RO is called from SUBROUTINE SOLVE. 

Thi3 subroutine i3 used to perform linear interpolation. SUBROUTINE 
LININT is called from SUBROUTINE RO. 

This subroutine is used to calculate rhe modified density coeffici- 
ents at the 5, H, and Q half points in the wing computational mesh. 
The modified density coefficients are computed by forming the approp- 
riate averages of the physical density coefficients calculated by 
SUBROUTINE RO. If Che flow is locally : upersonic, then upwinding is 
effected by forming a weighted average of the physical density at the 
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point in consideration and the density at a point in the upwind dir- 
ection along the appropriate curvilinear coordinate. Upwindlng io 
alunyo performed along the £ -wraparound coordinate. Upwinding in the 
H-opanwiae direction and in the £ -normal direction are uscr- 
defoatable options. SUBROUTINE ROCO ta called from SUBROUTINE SOLVE. 

RESID. Thia 3ubroutino la uoed to apply the full-potential equation residual 
operator to tho current velocity potential field to obtain the maos 
residuals at all points in the wing computational nc3h except for 
those points on the outer computational boundary. The mass residual 
at a given point io computed by forming the sum of the finite- 
difference derivatives of the mass fluxes in the q, n, and £ coordi- 
nate directions. Central difference expressions are used at interior 
field points while special procedures are used at the body surface. 
SUBROUTINE RESID employs the modified density coefficients calculated 
by SUBROUTINE ..0C0. SUBROUTINE RESID is called from SUBROUTINE 
SOLVE. 

OUTPUT . This subroutine is used to print the wing solution convergence 

history and to print the surface pressure coefficient distribution 
for each wing section. If desired, the density and Mach number con- 
tours may also be output. SUBROUTINES FORCE and TFORCE are called to 
compute the sectional and total force and moment coefficients, 
respectively. SUBROUTINE OUTPUT is called from SUBROUTINE SOLVE. 

GSTMET . This subroutine is used to retrieve the metric parameters from disk 
storage by in turn calling SUBROUTINE READMS . SUBROUTINE GETMET is 
called from SUBROUTINE SOLVE. 

FORCE . This subroutine is used to compute the wing sectional lift, 

wave-drag, and moment coefficients by integrating the surface 
pressure distribution for a given spanwise station. SUBROUTINE FORCE 
is called from SUBROUTINE OUTPUT. 

TFORCE . This subroutine is U3ed to compute the wing total lift and wave-drag 
coefficients. The moment coefficient at each spanwise station is 
also computed using the root quarter chord station as the datum about 
which Qonent3 are computed. For this subroutine to execute properly, 
the relations JSKIP-1 and JCPEND>NJIM must be satisfied. SUBROUTINE 
TFORCE is called from SUBROUTINE OUTPUT. 

READEC . This subroutine is used to transfer data from extended core storage 
to central memory (replaces CDC SUBROUTINE READEC). 

WRITEC . This subroutine is used to transfer data from central memory to ex- 
tended core storage (replaces CDC SUBROUTINE WRITEC) . 

READMS . This subroutine is used to transfer data from a mass storage device 
to central memory (replaces CDC SUBROUTINE READMS). 

NRITMS . This subroutine is used to transfer dat3 from central memory to a 
mass storage device (replaces CDC SUBROUTINE WRITMS). 

SAVE . This subroutine is used to save the flow field solucion on disk file 
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TAPE 12 for raetarting tha program. 

TRIB» Thio subroutine la uoed to solve a systea of simultaneous linear 

equations irtth a tridiagonal coefficient natrix and fixed boundary 
conditiono. Thio routine is uoed in the AF2 iteration algorithm to 
obtain a potential function intermediate correction. SUBROUTINE TRIB 
is called from SUBROUTINES HS0LV3 and SOLVE. 

TRIP. This subroutine is used to solve a systea of simultaneous linear 

equations with a tridiagonal coefficient matrix and periodic boundary 
conditions. Thi3 routine is used in the AF2 iteration algorithm to 
obtain a potential function intermediate correction. SUBROUTINE TRIP 
is called from SUBROUTINES NSOIVE and SOLVE. 


SECANT , This subroutine is used to determine the nacelle local internal flow 
Mach number using one-dimensional gas dynamic formulae which express 
mass continuity for the captured streamtube. A reference area and 
Mach number are supplied, and from the local flow area the local Mach 
number is determined. A secant iteration scheme is used to numeri- 
cally determine the zero of the governing implicit relation. SUB- 
ROUTINE SECANT is called from SUBROUTINE NINIT and determines the 
compressor face Mach number from the user-supplied capture ratio. It 
also determines the local Mach number from the compressor face Mach 
number and the ratio between the local and compressor face areas. 

METRIC . This subroutine is used to compute the metric parameters for all 

points in the nacelle component grid. The metric coefficients define 
the three-dimensional grid mapping from physical space to computat- 
ional space. The metric coefficients are computed using second-order 
accurate finite-difference formulae. The physical grid point locat- 
ions are entered through a binary read of a disk file. SUBROUTINE 
METRIC is called from SUBROUTINE NINIT. 

NRO . This subroutine is used to calculate the physical density coeffici- 
ents at the C-half points and the £-end points in the nacelle compu- 
tational me9h. From the current potential function field, the deriv- 
atives of the velocity potential in the £, n, and Q directions are 
computed using second-order accurate finite-difference formulae. The 
potential function derivatives are then used in a form of the isen- 
tropic energy equation to determine the density. At the body sur- 
face, the flow tangency condition is used in computing the density. 

At the internal outflow surface, the density is assumed known from 
the prescribed boundary conditions. At the outer computational 
boundary, the density is calculated using potential values inter- 
polated from the wing solution. SUBROUTINE NRO is called from SUB- 
ROUTINES NSOLVE and NOUTPUT. 

NROCO . This subroutine is used to calculate the modified density coeffici- 
ents at the 5,0, and C half points in the nacelle computational 
mesh. The modified density coefficients are computed by forming the 
appropriate averages of the physical density coefficients calculated 
by SUBROUTINE NRO. If the flow is locally supersonic, then upwinding 
is effected by forming a weighted average of the physical density 3t 
the point in consideration and the density at a point in the upwind 
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direction along the appropriate curvilinear coordinate. Upvinding io 
alvaya porformed along the ^-wraparound coordinate. Upvinding io not 
parforccd in the 0 - circumferential direction but can be in the C - 
radial direction. SUBROUTINE NROCO io called from SUBROUTINE NSOLVE 

This cubroutlne ia uoed to apply the full-potential equation residual 
operator to the current velocity potential field to obtain the mass 
reoidualo at all polnto in the nacelle computational nosh except for 
those points on the outer computational boundary. The maao reoidual 
at a given point io computed by forming the sum of the finite- 
difference derivatives of the maso fluxes in the £, n, and C coordi- 
nate directions. Central difference expressions are used at interior 
field points while special procedures are used at the body surface 
and at the internal outflow surface. SUBROUTINE NRESID employs the 
modified density coefficients calculated by SUBROUTINE NROCO. SUB- 
ROUTINE NRESID is called from SUBROUTINE NSOLVE. 

This subroutine is the main nacelle flow field integration control 
routine and is used to apply the AE2 approximate factorization scheme 
to determine the Inviscid nacelle flow solution. The AF2 scheme is 
applied in an iterative manner with the iteration sequence being 
terminated by either achieving convergence or by reaching a maximum 
permissible number of iterations. Per iteration, a three-step pro- 
cess is used to obtain the potential function correction. The first 
two steps calculate intermediate potential function corrections aris- 
ing from solving systems of factorized difference equations in the 
s-wraparound and n-circumferential directions. The third step calcu- 
lates the actual potential correction by solving a system of factor- 
ized difference equations for the ^-radial coordinate. Steps 1 and 2 
require sweeping the computational mesh from the body surface out- 
wards. The third step sweeps from the outer computational boundary 
towards the body surface. SUBROUTINE NSOLVE calls SUBROUTINES NRO, 
NROCO, and NRESID for calculation of the physical density 
coefficients, modified density coefficients, and mass residuals, 
respectively. SUBROUTINES TRIB and TRIP are called in the course of 
applying the AF2 scheme to solve systems of tridiagonal simultaneous 
linear equations with fixed and periodic boundary conditions, 
respectively. SUBROUTINE NOUTPUT is called to print the computed 
surface solution. SUBROUTINE NSOLVE is called by routine TWN. 

This subroutine is used to print the computed nacelle surface 
solution when convergence has been attained or when the maximum 
permissible number of iterations has been reached. SUBROUTINE NRO is 
applied at the body surface to calculate the physical density 
coefficients at mesh mid-points on the surface. From the computed 
density distribution, the surface Mach number and pressure 
coefficient distributions are determined. SUBROUTINE NOUTPUT is 
called from SUBROUTINE NSOLVE. 

This subroutine is used to transfer property data between central 
memory ar.d a. mass storage device. 
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IirFOT PARAMETERS AND OUTPUT INTERPRETATION 


1. INTRODUCTION 

In this section, cho input parameters are defined for the GRGSN3 and 
NGRIDA grid generation programs, the LK grid interface program, and the TWN 
wing/nacelie flow solution progrea. Where applicable, both the default value 
and the typical value of each input parameter are given. Discussions on 
interpretation of the output, PARAMETER 3tateaent specification, and the file 
usage for each program are also provided. 

2. GRGEN3 WING GRID GENERATION PROGRAM INPUT PARAMETERS 

The input data required for execution of the GRGEN3 wing grid generation 
program are entered by namelist input and, in some cases, formatted read 
statements. In all cases, the two NAMELISTS GRID3D and GRIDIN are entered. 
Most of the input parameters have assigned default values. 

2.1 NAMELIST GRID3D 

The input parameters entered in NAMELIST GRID3D specify the number of 
mesh points in the wing grid, control the wing grid mesh spacing, and Identify 
the location of the engine nacelle. 

NI A positive integer denoting the number of q -wraparound stations in 

the wing computational grid, a3 shown in Figure 25. The Sj-statlons 
are numbered in a counterclockwise direction starting at the upper 
surface trailing edge. The specified value of NI must be odd. A 
default value of 101 is specified for NI. 

NJ A positive integer denoting the number- of n~spanwlse stations in the 

wing computational grid, as shown in Figure 25. The n-stations are 
numbered starting at the wing root and increasing in the outward 
direction. A default value of 26 Is specified for NJ. 

NK A positive integer denoting the number of C-radial stations in the 

wing computational grid, as shown in Figure 25. The C-station3 are 
numbered starting at the outer boundary and increasing in the Inward 
direction. A default value of 18 is specified for NK. 

NJT A positive integer denoting the index of the first n-station that is 

off the wing tip. A default value of 19 is specified for NJT. 

NJC A positive integer denoting the index of the H-station that contains 

the nacelle plane of symmetry. A default value of 9 is specified for 
NJC. 

SFACT A real variable denoting the stretching factor used in controlling 

the n~spanwise grid spacing beyond the wing tip. A default value of 
1.2 is specified for SFACT. 
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Figure 25. Topology for wing cocponent grid. 
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SMOOTH A real variable denoting the factor used in smoothing the grid span- 
wise spacing near the wing tip. A default value of 1.0 is specified 
for SMOOTH. 

SVPANG A real variable denoting the wing sweep angle, in degrees, along the 
spanwise direction. A default value of 0.0 is specified for SWPANG. 

TWANG A real variable denoting the wing twist angle, in degrees, along the 
spanwise direction. A default value of 0.0 is specified for TWANG. 

CEL A real variable denoting the distance between the nacelle plane of 

symmetry and the wing root station. A default value of 3.0 is speci- 
fied for CEL. 

TIP A real variable denoting the spanwise location of the wing tip (the 

tip is assumed to lie halfway between the span stations NJT and 
NJT-l). A default value of 6.0 is specified for TIP. 

DYR00T A real variable denoting the spanwise distance betveen the wing 

second span station and the root station. If DYR00T * 0.0 is enter- 
ed, then DYRC0T is internally reset to (TIP/(NJT-1 . 5 ) ) . A default 
value of 0.0 is specified for DYR00T. 

DYTIP A real variable denoting the spanwise wing grid spacing near the wing 
tip. If DYTIP « 0.0 is entered, then DYTIP is internally reset to 
(TIP/(NJT-1.5)). A default value of 0.0 is specified for DYTIP. 

DYN’CL A real variable which controls the wing grid spanwise spacing near 
the nacelle plane of symmetry. A default value of 0.1 is specified 
for DYNCL. 

0.2 NAMELIST GRIDIN 

The irput parameters entered in NAMELIST GRIDIN specify the wing geometry 

and control the relaxation solution algorithm used to generate the wing grid. 

RADMAX A positive real variable denoting the radius of the circular outer 
boundary of the wing computational grid. A default value of *>.0 is 
specified for RADHAX. 

N3TATN An integer denoting whether one or more airfoil sections are to be 
used in defining the wing geometry. If NSTATN « 0, one airfoil 
section will be used. If NSTATN > 0, more than one airfoil section 
will be used. A default value of 0 is specified for NSTATN. 

I0PT An integer denoting the wing section geometry option. I0PT specifies 
the following geometry options: 

IOPT Option 

1 NACA OO.tX type airfoil 

2 circular arc airfoil (parabolic lead- 
ing edge blunting) 

3 circular cylinder 
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read In arbitrary airfoil geonatry 
circular arc airfoil (square root 
leading edge blunting) 

ONERA M6 airfoil 
KORN airfoil 


A default value of 1 ia specified for IOPT. 

TMAX A positive real variable denoting the airfoil maximum thickness as a 
fraction of the chord length. TMAX oust be entered only if IOPT * l, 
2, or 5. A default value of 0.12 is specified for TMAX. 

XC A real variable denoting the leading edge bluntness parameter. XC 

oust be entered only if IOPT ■ 2 or 5, 3nd must be in the range of 
0.0 < XC <0.4. A default value of 0.2 is specified for XC. 

XCN A real variable denoting the x-coordinate , in tons of fraction of 

the chord, of the point to which the origin is shifted. A default 
value of 0.0 is specified for XCN. 

3INX A real variable denoting the stretching parameter which controls the 
grid point distribution on the wing surface. 3 INN’ controls the sur- 
face grid point distribution as follows: 

3IXN Option 

0.0 the grid points -ire distributed evenly 

in terms of the arc length 
-1.0 the grid point distribution is 

established from a data statement which 
was developed by using a circle plane 
conformal mapping around a NACA 0012 
airfoil 

> 1.0 the grid point distribution is estab- 

lished from an exponential clustering 
formula with a BIN'N value closer to 
1.0 cluster ng more points at the lead- 
ing and trailing edges of the wing 
sections 

A value of 1.05 Is suggested for BIN’N. A default value of -1.0 is 
specified for BINN. 


ERCRI0 A real variable denoting the convergence tolerence used in the wing 
grid generation solution algorithm. A default value of 0.001 is 
specified for ERGRID. 

MAXIT An integer denoting the maximum number of iterations permissible to 
attain convergence in the wing grid generation algorithm. A default 
value of 200 is specified for MAXIT. 

AXCRI0 A real variable denoting the largest value in the 3lpha acceleration 
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parameter sequence used in the wing grid generation relaxation algo- 
rithm. A default valua of 0.01 is specified for AH GRID. 


* 

• ALGRID 

A real variablcj denoting the smallest value in the alpha acceleration 
parameter sequence used in the wing grid generation relaxation algo- 
rithm. A default value of 0.00006 is specified for ALGRID. 


- 

KGRID 

An integer denoting the starting element in the alpha sequence used 
for the wing grid generation. A default value of 0 is specified for 
KGRID. 



MGR ID 

An integer denoting the number of elements in the alpha sequence used 
for the wing grid generation. A default value of 3 Is specified for 
MGR ID. 



•DM EG 

A positive real variable denoting the relaxation parameter used in 
the wing grid generation algorithm. A default value of 2.0 Is speci- 
fied for DMEG. 



N0UT5 

An integer which controls the amount of grid output information. If 
N0UT5-1 , the grid point coordinates are printed. If N0UT5-0, the 
3 rid point coordinates are not printed. A default value of 0 Is 
specified for N0UT5. 



N'OUTo 

An integer denoting the increment in iteration count at which the 
grid point coordinates are to be printed during the course of the 
relaxation solution. A default value of 201 is specified for N0UT6. 
(If convergence is attained before N0UT6 iterations, the N0UT5 print 
option is invoked.) 

■ 


IINCR 

An integer which controls the grid solution output frequency in the E - 
direction. The wing grid coordinates will be printed with an IINCR 
increment. A default value of 3 is specified for IINCR. 

» 


JI.N'CR 

An integer which controls the grid solution output frequency in the n - 
direction. The wing grid coordinates will be printed with a JINCR 
increment. A default value of 6 is specified for JINCR. 



2.3 IN’ PUT DATA FOR WING MULTIPLE STATION DEFINITION OPTION 



The following input parameters are entered by formatted read statements 
only for the IOPT-4 wing multiple station definition option. No default 
values are specified for the following parameters. 



CARD L 

[A-Format 1 : 


m 

CHAR 

Title card (any alphanumeric characters). 

, 


CARD 2 

[Format (7F10.6)] : 


• 

Y PLANE 

A real variable denoting the local wing section span position (y 
coordinate) normalized by the root chord. 



CHORD 

A real variable denoting the local wing section chord normalized by 
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the root chord 
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THICK A real variable denoting the local wing section thickness in terns of 
fraction of the root chord (if the wing has the sane section but with 
varying thickness). 

ALPHA A real variable denoting the local wing section twist angle, in 
degrees. 

ALAMB A real variable denoting the local wing section sweep angle, in de- 
grees. 

XCEN A real variable denoting the x-coordinate , in tents of the fraction 
of the chord, to which the origin is shifted. 

RADIUS A real variable denoting the radius of the circular outer boundary. 

For highly tapered wings, the values of radius should be decreased in 
the spanwise direction. 

CARD 3 (Format (15) J: 

NPTS An integer denoting the nunber of points in the wraparound direction 
used to define the local wing section. The points are numbered from 
the trailing edge counterclockwise and do not include the redundant 
trailing edge point. 

CARD 4 [Format (2F10.5)): 

XB A real variable denoting the wing section chordwise coordinate, in 

terms of fraction of the local chord. 

Y3 A real variable denoting the wing section normal coordinate, in terms 

of fraction of the local chord. 

Note there are NPTS CARD 4 cards for each wing defining station. CARDS 

1-4 must be supplied for each of the NSTATN defining stations. 

2.4 PARAMETER STATEMENT SPECIFICATION 

The GRCEN3 wing grid generation program uses variable array dimensions. 

The respective array sizes are fixed by specification of the following dimen- 
sion parameters on the PARAMETER statement: 

Parameter Allowed Values 

IZ2 _> NI 

122 „ IZZ+2 

JZ > NJ 

KZ > NX 

The array dimensions will be fixed at the time of compilation. 





55 



2.5 


FILE USAGE 


\ " 


ORJG&’AL r.V: ...” ;.*> 
OF POOR quality 


The following flloa are used by Che GRGEN3 wing grid generation code. 
File No. Usage 


TAPE 5 
TAPE 6 
TAPES 9,10 


Input file 
printed output flic 
output flle3 containing the 
wing grid coordinates 


3. GRGEN3 PROGRAM OUTPUT INTERPRETATION 


The initial portion of the computer printout for the GRGEN3 wing grid 
generation program consists of the supplied input parameters. The remaining 
output parameters are defined below. 


J 

Y 

I 


X3.YB 


S2 

SX 


index of the spanwise wing sections 
final spanwise coordinate distribution 
index of points in the wraparound 
direction (starting from the upper 
trailing edge point) 

wing section coordinates (the coordinates 
are printed twice, first the input 
coordinates are printed and then the final 
coordinates after interpolation and 
clustering are printed) 
arc length normalized by section chord 
first difference of the arc length 


The following parameters are associated with the convergence history for the 
ADI numerical grid generation scheme. 


N 

CXMAX 

I.J 

CYMAX 

I.J 

RXMAX 

I.J 

RYMAX 

I.J 

ALPHA 


iteration number in the numerical grid 
generation scheme 

maximum correction (chordwise direction) 
position of the maximum correction CXMAX 
maximum correction (normal direction) 
position of the maximum correction CYMAX 
maximum residual (chordwise direction) 
position of the maximum residual RXMAX 
maximum residual (normal direction) 
position of the maximum residual RYMAX 
the acceleration parameter used 
for the ADI scheme 


4. NGRIDA NACELLE GRID GENERATION PROGRAM INPUT PARAMETERS 


The input data required for execution of the NGRIDA nacelle grid gener- 
ation computer program are entered by namelist input. In all cases, the four 
NAMELISTS LIST1, LIST2, LIST3, and LIST4 are entered. All of the input para- 
meters have assigned default values. A suitable grid can often be generated 
by retaining many of the input parameters at their default values. 
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Tha Input parameters specified In NAMELIST LIST1 specify the nunber of 
mesh points In the nccalle grid, outer boundary shape, grid point distribution 
stretching factors, and print and plot options. 

The following three parameters specify the number of grid points. 

IMAX A positive integer denoting the nunber of 5-wraparound stations In 

the nacelle computational grid, as shown In Figure 26. The specified 
value of IMAX must be less than or equal to 140, and should be of the 
fora (3cH-l) If grid sequencing is to be used In the GRAPE algorithm, 
where a is an integer. A default value of 67 is specified for IMAX. 

JMAX A positive integer denoting the number of ircumferential stations 
In the nacelle computational grid, as shown in Figure 27. Stations 
J-l and J“JMAX are coincident. The remaining meridional stations are 
spaced in equal angular increments around the body. The specified 
value of JMAX must be odd. A default value of 13 is specified for 
JMAX. 

KMAX A positive integer denoting the number of T-radial stations in the 
nacelle computational grid, as shown in Figure 26. The specified 
value of KMAX must be less than or equal to 60, and should be of the 
form (3n+l) if grid sequencing is to be used in the GRAPE algorithm, 
where n is an integer. A default value of 13 is specified for KMAX. 

The next three parameters specify the nacelle grid outer computational 
boundary. 


XLEFT A real variable denoting the x-coordinate value of the left-hand side 
of the nacelle grid outer computational boundary, as shown in Figure 
23. XLEFT is specified in the units of the input geometry, and typi- 
cally is negative. A default value of -30.0 is specified for XLEFT. 

RADOUT A positive real variable denoting the Radius of the nacelle grid 
outer computational boundary, as shown in Figure 23. RADOUT is 
specified in the units of the input geometry. A default value of 
50.0 is specified for RADOUT. 

DELTA A positive real variable denoting the offset distance between the 

centerline computational boundary and the x-axls, as shown in Figure 
28. DELTA is specified in the units of the input geometry and must 
not be zero. A value of approximately i% of the compressor face 
radius is recommended for DELTA. A default value of 0.5 is specified 
for DELTA. 

The next four parameters control the grid point distribution on the 
nacelle surface. 


IOUTER A positive Integer denoting the nunber of grid points desired on the 
outer portion of the nacelle surface ranging from the hillte point to 
the external outflow surface, as shown in Figure 2S. A default value 
of 33 is specified for IOUTER. 


IIN'NER A positive integer denoting the number of grid points desired on the 
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inner portion of the nacelle surface ranging from the hilite point to 
the compressor face outflow surface, as shown in Figure 28. The 
specified values for I0UTER and IINNER oust satisfy the relation 

IMAX =* IINNER + IOUTER - 1 (29) 

A default value of 30 is specified for IINNER. 

ALPHAO A positive real variable denoting the stretching factor used in 

determining the grid point distribution on the nacelle external sur- 
face. Entering ALPHAO=l.O produces a point distribution with uniform 
arc length between the points. Entering ALPHAO > 1.0 clusters points 
near the nacelle hilite. A default value of l.l is specified for 
ALPHA"!. 

ILP'LAI \ positive real variable denoting the stretching factor used in 

determining the grid point distribution on the nacelle internal sur- 
face. Entering ALPHAI=l.O produces a point distribution with uniform 
trc length between the points. Entering ALPHA I * 1.0 clusters points 
near the nacelle hilite. A default value of 1.1 is specified for 
.iLPHAI. 

The next five parameters determine the outer boundary grid point distri- 
bution. Generally, the outer boundary points near both the external and 
internal outflow surfaces are found by equating the axial coordinate of a 
given boundary point to the axial coordinate of the corresponding body point 
'see Figure 29). The remaining boundary points are found by either using an 
angular distribution or an arc-length distribution. 

IAA A positive integer denoting the wraparound coordinate station up to 

which the points on the computational boundary near the centerline 
have cheir axial coordinates equated to the respective values of the 
corresponding points on the internal nacelle surface, as illustrated 
in Figure 29. IAA is specified in terms of the GRAPE algorithm 
point-ordering scheme in which the wraparound coordinate initiates at 
the compressor face and terminates at the external outflow surface 
(see Figure 29). A default value of 12 is specified for IAA. 

IDD A positive integer denoting the wraparound coordinate station after 

which the points on the outer computational boundary have their axial 
coordinates equated to the respective values of the corresponding 
points on the nacelle external surface, as shown in Figure 29. IDD 
is specified in terms of the GRAPE algorithm point-ordering scheme in 
which the wraparound coordinate initiates at the compressor face and 
terminates at the external outflow surface (see Figure 29). Note 
that IDD must satisfy the relation IAA < IDD < IMAX. A default value 
of 57 is specified for IDD. 

■OUTER An integer denoting whether the outer boundary grid points between 

points LAA and IDD are to he determined by in angular distribution or 
by an 3rc-length distribution. If K9l'TKR=0, the grid points are 
determined by an angular distribution with the grid points being 
spaced at equal angular increments around the nacelle hilite. If 
OLTER=i, the grid points are determined by an arc-length distri- 



RADIUS IG=IDD-l IG=IDD 



Figure 29. Definition of points IAA and IDD. 



'I 


CRiK'j.J.i!, r s 'X'7 
OF POOH QUALITY 

bution with respective stretching factors being specified by the 
user* A default value of 0 is specified for KOUTER. 

If K0UTER“1 io specified, then the following two parameters must be entered. 

ALPHBO A positive real variable denoting the stretching factor used in 

determining the nacelle outer boundary grid point distribution for 
the portion of the boundary between points C and IDO, as shown In 
Figure 29. Entering ALPHB0-1.0 produces a point distribution with 
uniform arc length between the points. Entering ALPHBO > t.O 
clusters points near point C. A default value of 1.0 is specified 
for ALPHBO. 

ALPHBI A positive real variable denoting the stretching factor used in 

determining the nacelle outer boundary grid point distribution for 
the portion of the outer boundary between points C and IAA, as shown 
in Figure 29. Entering ALPHBI*!. 0 produces a point distribution with 
uniform arc length between the points. Entering ALPHBI > 1.0 clus- 
ters points near point C. A default value of 1.0 is specified for 
ALPHBI. 

The next three parameters control the grid spacing adjacent to the body 
surface in the direction normal to the surface. Generally, the normal 
distance DS , as shown in Figure 30, is computed by the relation 

DS = (RCF - DELTA)/ (UMAX - 1) (30) 

where RCF is the compressor face radius. For isolated nacelle grids, it is 
generally preferable to adjust DS near the external outflow surface to make 
the radial point distribution uniform as shown in Figure 30. This point dis- 
tribution produces a more favorable cell aspect ratio at the outflow surface, 
and thereby enhances stability and convergence in the flow solution algorithm. 

DSFACT A positive real variable used as a multiplier for DS which is defined 
by equation (30). A default value of 1.0 is specified for DSFACT. 

KRAMP An integer denoting whether or not the DS distribution is to be 

altered to produce a uniform radial point distribution at the extern- 
al outflow surface. If KRAMP*!, the DS distribution is altered to 
produce the uniform radial point distribution. If KRAMP=0, the DS 
distribution is not altered, and the grid points will be packed 
towards the nacelle surface. If an isolated nacelle grid is being 
generated, then the KRAMP*L option is recommended. A default value 
of 1 is specified for KRAMP. 

IRAMP A positive integer denoting the wraparound station after which the DS 
distribution is altered to produce a uniform radi3l point distri- 
bution at the external outflow surface, as shown in Figure 30. IRAMP 
must be entered only if KRAMP » 1 . IRAMP is specified in terms of the 
GRAPE algorithm point-ordering scheme in which the wraparound coordi- 
nate initiates at the compressor face and terminates at the external 
outflow surface. Note that IRAMP must satisfy the relation IRAMP < 
IMAX. A default value of 57 is specified for IRAMP. 
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Figure 30. Nacelle grid normal mesh spacing and 
definition of point IRA11P. 
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The next three parameters specify the scaling and smoothing function 

parometcro. 

SCALE A positive real variable denoting the scaling multiplier used in 

obtaining the nacelle final grid. Entering SCALE»l.O produces a grid 
in totno of the original input units. Setting SCALE“l.O/RCF, where 
RCF ia the compressor face radius, produces a grid with a compressor 
face radius of 1.0. This option is recommended as it allows for more 
easily determining the optimum acceleration parameters in the flow 
solution code. A default value of 0. 11904762-1 .0/8.4 is specified 
for SCALE. 

KTIME An integer denoting the number of times smoothing polynomials are 

applied in determining the nacelle grid outer boundary axial coordi- 
nate distribution in the vicinity of point IDD, and the DS distri- 
bution in the vicinity of point IRAMP if the KRAI IP =1 option is speci- 
fied. If KTIME =0 is entered, no smoothing is performed. The default 
and recommended value for KTIME is 3. 

MDEL A positive integer denoting the number of points to the left of point 
IDD at which the smoothing is initiated for the nacelle grid outer 
boundary point axial locations if KTIME > 0 is specified. NDEL must 
be entered only if KTIME > 0 is specified. A default value of 5 is 
specified for MDEL. 

The next two parameters specify the print options. 

KDUMP An integer denoting whether or not detailed nacelle surface and outer 
boundary coordinate data are to be printed. If KDUMP”0, the data are 
not printed. If KDUMP’ , 1 , the printout is performed. A default value 
of 1 is specified for KDUMP. 

KPRINT A one-dimensional integer array consisting of JMAX elements. Each 
element of KPRINT specifies whether or not a meridional plane grid 
point distribution is to be printed. Specifying KPRINT(J)=*1 
( J»1 , . . . , JMAX) causes the grid point coordinates to be printed for 
the Jth meridional plane. If KPRINT(J) a 0 is entered, the printout is 
skipped for the Jth plane. A default value of KPRINT(l)“l is 
specified, while the remaining elements of KPRINT are specified as 0. 


The next two parameters specify the plot options. 


KPL0T An integer denoting whether or not the J»1 meridional plane grid is 
to be plotted. Specifying KPL0T»1 causes the grid to be plotted for 
the J»1 meridional plane. If KPLOT’O is entered, the plotting is not 
performed. A default value of KPL0T”! is specified. 


KDEV 


An integer denoting the plotting device. KDEV must be specified only 
if KPL0T-1 is entered. If KDEV=l, the plot device is the VERSATEC 
electrostatic plotter. If KDEV-2, the plot device is the HP ink 
plotter. A default value of 1 is specified for KDEV. (Note that the 
plot device specification can be altered by modifying SUBROUTINE 


CONVERT). 
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Tha next parcsattr specifies chs taps number (disk file) on which the 
grid coordinate data are written. 

ITAPS An integer denoting tho tape (file) number on which the grid 

coordinate data ore written. ITAPE can ba entered as either 10 or 
14. For isolated nacelle computations, ITAPE should bo retained at 
its default vai of 10. For use of the nacelle grid in the TWN 
code, ITAPE ■ i4 * aould be entered. 

4,2 NAMELIST LIST2 

The input paraneters entered In NAMELIST LIST2 specify the nacelle/ inlet 
geometry which is described by tabular input. A default ax i symmetric 
nacelle/ inlet geometry has been loaded onto DATA statements within the 
program. The default geometry data are for the Lockheed-Oeorgia GELACi 
Jxisymmetric nacelle/ inlet configuration. 

The following three paraneters specify the internal surface contour for 
axisynmetric nacelles. 

\I A positive integer denoting the number of tabular data points used in 

defining the internal surface contour, as shown in Figure 31. A 
default value of 60 is specified for h'l. 

XIX A one-dimensional real variable array consisting of NI elements. 

XIN(I) (I»l,...,NI) denotes the axial coordinate of the Ith point 
used in defining the internal surface contour, as shown in Figure 31. 
Successive elements of XIN must be monotonically increasing. XIN(l) 
denotes the axial station of the nacelle hilite, and XIN(NI) denotes 
the compressor face axial station. Default values for XIN(I) 

(1*1, . .. ,60) are specified. 

XIN A one-dimensional real variable array consisting of NI elements. 

XIN(I) (1*1, ...,NI) denotes the radius of the Ith point used in 
defining the internal surface contour, as shown in Figure 31. RlN(i) 
denotes the nacelle hilite radius, and R1N(N1) denotes the compressor 
face radius. Default values for RIN(I) (1*1,..., 60) are specified. 

The next three parameters specify the external surface contour for 
axisyametric nacelles. 

NO A positive integer denoting the number of tabular data points used in 

defining the external surface contour, as shown in Figure 31. A 
default value of 100 is specified for NO. 

X0UT A one-dimensional real variable array consisting of NO elements. 

X0UT(I) (1*1, ...,N0) denotes the axial coordinate of the 1th point 
used in defining che external surface contour, as shown in Figure 31. 
The successive elements of X0UT must be monotonically increasing. 
X0UT(l) denotes the nacelle hilite axial station, and XOUT(NO) 
denotes che axial station of the external outflow surface. The 
relation XDL'T ( 1 ) * XIN(l) must be satisfied. Default values for 
XOl'T(I) (1*1 , . . . , 100) are specified. 
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A ona-dimensional real variable array consisting of HO elements. 
R0UT(1) (1®1,...,K0) denocoa the radius of the Tth point used in 
defining the external surfaca contour, as shown in Figure 31. 

R0UT(1) danotea Cha nacelle hilite radius, and ROUT(NO) denotes the 
nacelle radius at the external outflow surface. The relation 
ROUT(l)=RIN(l) must ba satisfied. Default values for RO’JT(I) 

(1*1 , . . . , 100) are specified. 

4.3 NAMELIST LIST3 

The input parameters entered in NAMELIST LIST3 specify the convergence 
criteria, iteration Units, and print options used in the GRAPE algorithm. 
These parameters are given only cursory descriptions herein, but are more 
fully discussed in Reference 6. Generally, the program is executed by 
retaining all of the parameters in NAMELIST LIST3 at their default values. 

NORDA A cne-dinens ional integer array consisting of two elements which 

denote the convergence criteria for the coarse-mesh and the fine-nesh 
GRAPE solutions. NORDA(l) and N0RDA(2) denote the numbers of order 
of magnitude by which the maximum correction is to be reduced for the 
coarse-mesh and fine-mesh solutions, respectively. Default values of 
NCRDA(l)-6 and \’ORDA'2)-3 are specified. 

MAXITA A one-dimonsional integer array consisting of two elements which 

denote the iteration limits for the coarse-mesh and fine-nesh GRATE 
solutions. MAXITA(l) and MAXITA(2) denote the iteration limits for 
the coarse-mesh and fine-mesh solutions, respectively. The MAXITA 
array is used in conjunction with the NORDA array. Default values of 
MAXITA( l )"400 and MAXITA (2 )»200 are specified. 

JPRT An integer denoting the print option to be used in the GRAPE 

algorithm. If JPRT < 0, no printing is performed except for error 
warning messages. If JPRT"l, a detailed printout will be performed. 

A default value cf -1 is specified for JPRT. 

4.4 NAMELIST MST4 

The input parameters entered in NAMELIST LIS74 specify the grid control 
functions and convergence parameters used in the GRAPE algorithm. Again, only 
cursory descriptions are provided herein, with detailed definitions being 
available In Reference 6. Generally, the program is executed by retaining all 
of the parameters In NAMELIST LIST4 at their default values. 

AAAI Positive real variables which control the enforcement of grid 

BBBI orthogonality and normal mesh spacing at the nacelle surface. Small 

values (i.e., 0.2) cause these effects to be propagated far into the 
field, whereas larger values (i.e., 0.6) cause these effects to 
diminish more rapidly, A default value of 0.6 is specified for both 
AAAI and 3BBI. 

CCCI Positive real variables which control the enforcement of grid 

DDDI orthogonality and normal mesh spacing at the outer boundary. Small 

values (i.e., 0.2) cause these effects to be propagated far into the 
field, whereas larger values (I.e., 0.6) cause these effects to 
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dininish core rapidly* A default value of 0.6 is apecified for both 
CCCI and DDDI. 

OKEGA A poaitivo real variable uaed in the GRAPE 

succcaoive-line-over-relaxation schema to obtain the grid point 
diotribution. OKEGA nuat be in the range of 0.0 < OMEGA < 2.0. A 
default value of 1*3 is specified for OMEGA. 

OMEGP Real variables uaed aa relaxation paraaetera in obtaining the body 
OMEGQ surface forcing functions that are uaed in the Poisson equations. 

The effects of controlling grid spacing and orthogonality at the 
nacelle surface car. be eliainated by entering OMEGP » OMEGQ ■ 0.0. 
OMEGP and OMEGQ must be in the range of 0.0 to 2.0. A default value 
of 0.1 is specified for both OMEGP and OMEGQ. 

OMEGR Real variables used as relaxation paraaeters in obtaining the outer 
0MEGS boundary forcing functions that are used in the Poisson equations. 

The effects of controlling grid spacing and orthogonality at the 
outer boundary can be eliminated by entering OMEGR » OMEGS ■ 0.0. 
OMEGR and OMEGS oust be in the range of 0.0 to 2.0. A default value 
of 0.0 is specified for both OMEGR and OMEGS. 

PLIM Real variables used as limitation factors in solving for the Poisson 
QLIM equation forcing functions. The range for each of these parameters 
RLIM is 0.0 to 100. 0. A default value of 0.5 is specified for PLIM, QLIM, 
SLIM RLIM, and SLIM. 

DSOBI A positive real variable denoting the normal distance to be imposed 

between the nacelle radial stations KMAX and KMAX-l . DSOBI has units 
of the input geometry. A default value of 0.2 is specified for 
DSOBI. 

4.5 PARAMETER STATEMENT SPECIFICATION 

The NGRIDA nacelle grid generation program uses variable array 
dimensions. The respective array sizes are fixed by specification of the 
following dimension parameters on the PARAMETER statement: 

Parameter Allowed Values 


NX 

> 

IM.JC 

NY 


l 

NYG 

> 

JMAX 

NZ 

> 

KMAX 

NX0 

> 

NO 

NXI 

*> 

NI 

NX0N 

> 

I0UTER 

NXIN 


I INNER 

NX02 


NO*NO 

NXI 2 

~> 

NI*NI 


The array dimensions will be fixed at the time of compilation. 
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The following files are used by the NGRIDA nacelle grid generation code. 


Flip Mo. 

TAPE 5 
TAPE 6 
TAPE 10 
or 

TAPE 14 
TAPES 1,4,7 


5. NGRIDA PROGRAM OUTPUT INTERPRETATION 


Usage 

input file 

printed output file 

output files containing the nacelle 

grid coordinates 

scratch files (should be deleted 
after computation) 


The initial portion of the computer printout for the NGRIDA program 
consists of the NAMELIST input data. Then if the KDUMP-L option is specified, 
detailed nacelle surface and outer boundary coordinate data are listed. After 
this, the GRAPE algorithm output is printed. The finished grid point 
coordinates are then printed far selected meridional planes as specified by 
the KPRINT array. Finally, certain index parameters are printed. 

The output parameters for the GRAPE algorithm are discussed in Reference 
6. The remaining output parameters are defined below. 


I £ -wraparound point index 

J n - civcumferential point index 

K r-radial point index 

X x-coordinate of point in finished 

grid, (scaled units) 

V y-coordinate of point in finished 

grid, (scaled units) * 

Z z-coordinate of point in finished 

grid, (scaled units) * 

XOUT x-coordinate of input point on external 

surface, (original units) 

ROUT radius of input point' or. external surface, 

(original units) 

XIN x-coordinate of input point on internal 

surface, (original units) 

RIN radius of input point on internal surface, 

(original units) 

ARC-LENGTH arc length measured from hliite point, 

(original units) 

XBO x-coordinate of redistributed point on 

external surface, (original units) 

RBO radius of redistributed point on external 

surface, (original units) 

XBI x-coordinate of redistributed point on 

internal surface, (original units) 

RBI radius of redistributed point on internal 

surface, (original units) 

X3URF surface point x-coordinate used as input 

fer GRAPE code, (original units) 

(Note - the y and z axes are along the wingspan for the wing and 
nacelle algorithms, respectively; the z and y smes are normal to 
the wing for the wing and nacelle algorithms, respectively) 
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RSURF surface point radius used as input for 

GRAPE code, (original units) 

XBOUND outer boundary point x-coordinate used as 

input for GRAPE code, (original units) 

RBOUND outer boundary point radiu3 used as input 

for GRAPE code, (original units) 

6. LK GRID INTERFACE PROGRAM INPUT PARAMETERS 

The input data required for execution of the LK grid interface program 
are entered by namelist input. In all case3, the single NAMELIST NACEL must 
be entered. 

6.1 NAMELIST NACEL 

The input parameters entered in NAMELIST NACEL specify the location of 
the nacelle with respect to the wing. 

XR A real variable denoting the distance, in terms of the fraction of 

the local wing section chord, between the wing leading edge and the 
nacelle hilite, as shown in Figure 32. A default value of 0.7 is 
specified for XR. 

ZT A real variable denoting the distance, in terms of the fraction of 

the local wing section chord, between the top of the nacelle and the 
point on the wing leading edge, as shown in Figure 32. A default 
value of 0.4 is specified for ZT. 

SCALE A real variable denoting the multiplying factor used to scale the 

nacelle and its grid. If the nacelle and wing reference lengths are 
the same, then SCALE^l.O. A default value of 0.22 is specified for 
SCALE. 

6.2 PARAMETER STATEMENT SPECIFICATION 


The LK grid interface program uses variable array dimensions. The 
respective array sizes are fixed by specification of the following dimension 


parameters on the PARAMETER statement. 

Parameter 


IZZ 2 

IZ2 

JZ _> 

KZ > 

IZJZ 
TSSZZZ 

IZT > 

IZB ^ 

NX > 


Allowed Values 
(see documentation for) 

NI (GRGEN3) 

IZZ-:-2 

NJ (CRGEN3) 

NK (GRGEN3 ) 

IZZ*JZ 

IZZ*JZ*KZ*6 

estimated number of points in the 
grid overlap region 
estimated number of points in the 
exhaust jet region at the nacelle 
aft station 

IMAX (NCR IDA) 
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Figure 32. Definition of nacelle position parameters XR 
and ZT. 



V « 


NY > 

NZ 2 

NXP 

TSS 

The array dimensions will be fixed at the 
6.3 FILE USAGE 

The following files are used by 

File No. 

TAPE 5 
TAPE 6 
TAPES 9,10 


TAPE 14 


TAPES 19,44 


7. LK PROGRAM OUTPUT INTERPRETATION 
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JMAX (NGRIDA) 
KMAX (NGRIDA) 
NX+l 

NX*NY*NZ*7 
tine of coapilation. 


LK grid interface program. 

Usage 

input file 

printed output file 

input files from the GRGEN3 wing 

grid generation code containing the 

wing grid coordinates 

input file from the NGRIDA nacelle 

grid generation code containing the 

nacelle grid coordinates 
output files from the LK grid 
interface program to be used as 

input to the TWN code 


The initial portion of the computer printout for the LK grid interface 
program consists of the NAMELIST input data. The remaining output parameters 
are defined below. 


XR 

XL 

ZT 

ZB 

XWOR 

YWOR 

ZWOR 

J1,J2 

NJC 

K1,K2 

KJ1.KJ2 

N 

NS 


the horizontal distance from the nacelle hilite to the wing leading 
edge (expressed in the wing coordinate system) 

Che horizontal distance from the nacelle aft station to the wing 

leading edge (expressed in the wing coordinate system) 

the vertical distance from the nacelle hilite (top meridional plane) 

to the wing leading edge (expressed in the wing coordinate system) 

the vertical distance from the nacelle hilite (bottom meridional 

plane) to the wing leading edge (expressed in the wing coordinate 

system) 

the x coordinate of the nacelle origin (expressed in the wing 
coordinate system) 

the y coordinate of the nacelle origin (expressed in the wing 
coordinate system) 

the z coordinate of the nacelle origin (expressed in the wing 
coordinate system) 

the J indices (spanwise direction) of the overlap region Inner 
boundary 

the J Index of che wing span station that contains the nacelle plane 
of symmetry 

the K indices (normal direction) of the overlap region inner boundary 
the K indices of the exhaust region at the aft station of the nacelle 
number of points in the overlap region (J’. > J < J2) 
number of points in the overlap region (J=«J1 or J»J2) 
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NN number of points In the exhaust jet region at the aft station of the 

nacollo 

8. TON FLOW SOLUTION PROGRAM INPUT PARAMETERS 

The input data required for execution of the TWN flow solution program 

are entered by namelist input. In all cases, the two NAMELISTS FLOWIN and 

FLOWC are entered. Most of the input parameters have assigned default values. 

8.1 NAMELIST FLOWIN 

The input parameters specified in NAMELIST FLOWIN control the computation 

of the wing component flow field. 

MODE An integer which controls the program mode of operation for the wing 
component flow field calculation. If M0DE»l, the mode of operation 
is for a wing/nacelle solution. If MODE»0, the mode of operation is 
for an isolated wing solution. A default value of 1 is specified for 
MODE. 

KCX An integer denoting the number of solution cycles between the wing 

and nacelle algorithms that are to be performed for wing/nacelle flow 
field computations. KCX =*l must be entered if M0DE»0 is specified. 

A default value of 4 is specified for KCX. 

FMACH A real variable denoting the free-stream Mach number. A default 
value of 0.75 is specified for FMACH. 

ALPHAW A real variable denoting the angle of attack, in degrees. A default 
value of 0.0 is specified for ALPHAW. 

ARAT A real variable denoting the wing aspect ratio based on a simple 

trapezoidal planfona. If ARAT=*0.0 is entered, the aspect ratio is 
computed internally. A default value of 0.0 is specified for ARAT. 

TIP A real variable denoting the span location of the wing tip. If 

TIP=0.0 is entered, then TIP is internally reset to 0.5* 
(Y(NJT)+Y(NJT-l)). A defaulc value of 0.0 is specified for TIP. 

ERROR A real variable denoting the convergence tolerance used for 

terminating the wing flow field calculation. Convergence is attained 
If the following relation is satisfied 

KaxUlCx! 1 CRH (3D 

where R^ ax is the maximum residual on the first iteration and is 

the maximum residual on the nth iteration. A default value of Q.uOl 
is specified for ERROR. 

MAXIT A positive integer denoting the maximum number of iterations per 

cycle permissible in calculating the wing component flow field. A 
default value of 40 is specified for MAXIT. 

The next four parameters specify the acceleration constants used in the 
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AF2 Iteration scheme for the wing algorithm. One aethod has been incorporated 
into tha progrca for calculating the acceleration parameter sequence. The 
method employs the sequence 

K-l 

o^' - Og (c^ / Ojj) M-l (K-1,2,3,...,M) (32) 

where a. is the lower limit fora , Og is the upper Unit for ct, k is the 
sequence element number, and M is the number of elements in the sequence. 

ALPHAL A positive real variable denoting the ct lower limit ct^ . 

A default value of 0.4 is specified for ALPHAL. 

ALPHAH A positive real variable denoting the ci upper limit ^g . 

A default value of 4.0 is specified for ALPHAH. 

M A positive integer denoting the number of elenents M of the sequence 

in equation (32). A default value of 3 is specified for >1. 

KKK An integer denoting the starting element of the sequence given in 

equation (32). KKK must satisfy the relation C MOCK M. A default 
value of 0 is specified for KKK. 

The following parameter specifies the standard relaxation factor used in 
the AF2 iteration scheme for the wing algorithm. 

OMEGA A positive real variable denoting the wing algorithm relaxation 

factor. For stability, OMEGA must be in the range 0.0 < OMEGA < 2.0. 
The default and recommended value for 0MECA is 1.8. 

The next three parameters control upwinding of the density and the amount 
of artificial viscosity and damping used in stabilizing the computations in 
regions of supersonic flow in the wing flow field solution. Mote that 
upwinding is always applied in the ^-wraparound coordinate direction. 

MDIF An integer which controls the upwinding of the density coefficients 
in the wing algorithm for regions of supersonic flow. MDIF controls 
the upwinding options as follows: 

MDIF Option 

0 upwinding in the wraparound direction 

1 upwinding in the wraparound and spanwise 
directions 

2 upwinding in the wraparound and normal 
directions 

3 upwinding in the wraparound, spanwise, and 
normal directions 

In most cases, upwinding is not required in all three directions. A 
default value of 1 is specified for MDIF. 

COM A positive real variable denoting the artificial viscosity 

coefficient used in performing upwinding of the density in regions of 
supersonic flow. CON generally takes the values 1.0 ^_C0N 2.0, 
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with higher values of CON producing greater upwinding. High values 
of CON can decrease shock wave recolution. A default value of 1.2 is 
specified for CON. 

A positive real variable denoting the tioelike dissipation 
coefficient which is used in the AF2 algoritha to produce diagonal 
dominance in the ^-difference equations and thereby enhance stability 
in regions of supersonic flow. Increasing BXI increases the amount 
of dissipation. A default value of 0.1 is specified for BXI. 

A real variable used as a circulation relaxation parameter in the 
wing flow field solution algorithm. Low values of RGAM retard the 
build up of circulation whereas higher values have the opposite 
effect. Too high of values of RGAM can cause divergent oscillations 
in the iterative solution process. Best results are generally 
achieved using a high RGAM value (1.3 ~ 1.5) for the first several 
Iterations and then reducing RGAM to a lower value (0.7 ~ 0.9) for 
the remaining iterations. A default value of 0.7 is specified for 
RGAM. 

An integer which denotes the interval in iteration number at which 
the wing solution is printed. If INTOUT-O, the printing of the 
solution is performed only at the termination of the calculation. A 
default value of 0 is specified for INTOUT. 

An integer which controls the amount of printout for the wing flow 
field solution algoritha. IOUT controls the printout as follows: 

IOUT Print/Plot Option 

0 density, Mach number and pressure 
coefficient plots 

1 density and Mach nunber plots 

2 density and pressure coefficient plots 

3 density plot 

4 pressure coefficient plot 

5 Mach number and pressure coefficient plots 

6 Mach number plot 

7 skip all print routines 

A default value of 4 is specified for IOUT. 

An integer denoting the last H-spanwise wing station up to which the 
wing solution is printed. If JCPEND-0 is entered, then JCPEND is 
internally reset to NJT-1. A default value of 0 is specified for 
JCPEND. 

An integer denoting the increment in n-spanwise station nutnber at 
which the wing solution is output. A default value of l is specified 
for JSKIP. 

An integer denoting whether or not to store the computed surface 
pressure coefficient data on disk file TAPE 20. If MCP*0, the file 
write is performed. If MCP=0, the file write is not performed. A 
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MORIN An Integer employed In restarting the program from a preceding 

calculation. If KORIN^O, the potential field is read from disk 
file TAPE 11. If M0RIN“O, the file read is not performed. A default 
value of 0 is specified for MORIN. 

MOROUT An integer employed in storing the potential field for 

post-processing or a restart calculation. If MOROUT^O, the 
potential field is written onto disk file TAPE 12. If M0R0UT**0, the 
file write is not performed. A default value of 0 is specified for 
MOROUT. 

After NAMELIST FLOWIH is entered, a title card is then entered. 

8.2 NAMELIST FLOWC 

The input parameters specified in NAMELIST FLOWC control the computation 

of the nacelle component flow field. 

XMCF A real variable denoting the effective Mach number at the compressor 
face. XMCF must be entered only if the inlet capture ratio (CRATIO) 
is not entered. A default value of 0.3516 is specified for XMCF. 

CRATIO A real variable denoting the inlet capture ratio, which is defined by 

CRATIO = AINF/AHL (33) 

where AINF is the cross-sectional area of the capture streamtube far 
upstream of the inlet, and AHL is the effective hiltr.e area (see 
Figure 33). If CRATI040-0, then the compressor face Mach number 
XMCF will be internally computed from the supplied value of CRATIO. 

If CRATIO ■ 0.0, then XMCF must be entered directly. A default value 
of 0.0 is specified for CRATIO. 

KINIT An integer denoting the nacelle grid potential function 

initialization option. If KINIT * 0, the potential function 
initialization is performed using free-streaa velocity components for 
each point in the computational mesh. If KINIT =* 1, the potential 
function initialization is performed using free-stream velocity 
components for the external flow mesh points, and using velocity 
components calculated from one-dimensional gas dynamic formulae for 
the Internal flow mesh points. The local internal flow Mach number 
is calculated using the compressor face Mach number and the ratio 
between the local flow area and the compressor face area. Setting 
KINIT =* 1 generally significantly improves convergence speed. The 
default and recommended value for KINIT is 1. 

XMITLM A positive real variable used as a Mach number multiplier for the 
nacelle grid internal flow points if che KINIT = 1 initialization 
option is specified. XMULM must be entered only if KINIT = 1, and 
generally is in the range of 0.6 to 0.3. The default value for XMULM 
is 0.7. 

CRIT A real variable denoting the convergence criterion used for 

terminating the nacelle flow field calculation. Convergence is 
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attained if the following relation is satisfied 

| B S a *|/| 8 i«| <. CR IT (34) 

x r* 

where ia the maximum residual on the first iteration, and R^ aj; 

is the maximum residual on the nth iteration. A default value oV 

0.001 ia specified for CRIT. 

ITMAX A positive integer denoting the maxima nunoer of iterations per 

cycle permissible in calculating the nacelle component flow field. A 
default value of 50 is specified for ITKAX. 

The next four parameters specify the acceleration constants used in the 
AF2 iteration scheme for the nacelle algorithm. Two options nave been 
incorporated into the nacelle program for calculating the acce •. e ration 
parameter sequence. The first option employs the sequence 


K-l 

\ * 0t H (' t L / a ri^ M " L (K-l, 2, 3, . . .,M) 


(35) 


where cij is the lower limit for 2, is the upper limit for a, K is the 
sequence' element number, and M is the number of elements in the sequence. The 
second option employs the sequence 


*1 * \ 


(36) 


~i * cc 

K K-l 

in which case M equals 8. 


28 


(9-K) (K-2,3, . ..,7,8) 


KALPHA An integer denoting the a-sequence to be used in the computation. If 
KALPHA-1 , the sequence given by equation (35) is used. If KALPHA-2 , 
the sequence given by equation (36) is used. The default and 
recommended value for KALPHA is 1. 


ALPHAL A positive real variable denoting the a lower limit n 
A default value of 1.0 is specified for ALPHAL. 

ALPHAH A positive real variable denoting the a upper limit a 
A default value of 20.0 is specified for ALPHAH. 

MD A positive integer denoting the number of elinents M of the sequence 

in equation (35). A default value of 3 is specified for HD. HD must 
be entered as 3 if the KALPHA-2 option is specified. 

The following parameter specifies the standard relaxation factor used in 
the AF2 iteration scheme for the nacelle algorithm. 


0HEGA A positive teal variable denoting the nacelle algorithm relaxation 

factor. For stability, OMEGA must be in the range 0.0 < OMEGA < 2.0. 
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Tha default and recommended value for OMEGA is 1*0. 

The next three parameters control the upwinding of tha density end the 
amount of artificial viscosity and dissipation used in stabilising tha 
computations in regions of supersonic flow in tha naccllo flow field solution. 

< Note that upuinding is always applied in tha 5 -wraparound coordinate 
direction. 

NKVIS An integer denoting whether or not upvlnding in regions of supersonic 
flow is to be applied in tha ^-radial coordinate direction. If 
NKVIS'O, the upwinding is not applied. If NKVIS»1, upwinding of the 
density is applied in the ^-direction. A default value of 0 is 
specified for NKVIS. 

CFACT A positive real variable denoting the artificial viscosity 

coefficient used in performing upwinding of the density in regions of 
supersonic flow. CFACT generally takas the values 1.0 «■ CFACT <_2.0, 
with higher values of CFACT producing greater upwinding.' High values 
of CFACT can decrease shock wave resolution. A default value of 1.2 
is specified for CFACT. 

3XIS A positive real variable denoting the tiaelike dissipation 

coefficient which is used in the AF2 algorithm to produce diagonal 
dominance in the ^-difference equations and thereby enhance stability 
in regions of supersonic flow. Increasing BXIE increases the amount 
of dissipation. A default value of 0.1 is specified for BXIE. 

The next three parameters control the smoothing operations performed in 
the computations. 

KSMTH A positive integer denoting the 4-radial station index K up to which 
smoothing is performed for the computed density in the •"-wraparound 
coordinate direction in the vicinity of the ; TRANE £-s~atior» (see 
Figure 34). The density smoothing generally Improves convergence 
speed. A default value of 4 is specified fov KSMTH. 

KTIME An integer denoting the number of times a smoothing polynomial is 
applied in smoothing the density values In the E-directlon in the 
vicinity of the ITRAN2 ^-station for the 4-stations ranging from l to 
KSMTH. If KTIHE-0 is entered, no smoothing is performed. The 
density smoothing generally improves convergence speed. A default 
value r 3 is specified for KTIME. 

NDEL A positive integer denoting the number of ^-stations to the left and 
to the right of point ITRAN2 at which the density is smoothed. 
Smoothing is performed from point ( ITAAN2-NDEL) to point 
(ITRAN2+-NDEL) for the 4 -stations ranging from 2 to KSMTH, and from 
point (ITRAN2+1) to point (ITRAN2+2*NI>SU for the 1st 4-station. A 
default value of 3 is specified for NEEL. 

The next two parameters specify the potential function under-relaxation 
factors. 


CORFAT A one-dimensional real variable array consisting of KMAX elements. 
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CORFAT(K) (K«1,...,KKAX) 9pcclfla3 tho potential function 
und or- relaxation factor for the Kth ^-station. This undcr-rolaxation 
to nacoaoary to oaintain stability and enhance convergence for mesh 
points near the Inlet centerline. Under-relaxation of the potential 
function is performed only for the C-atatlons which are in the range 
ITRAH3 <I< EL\X (see Figure 34). The eleaents of CORFAT generally 
range from "a low value at K * 1 to a value of 1.0 (no 
under-rclaxatlon) as K approaches KKAX. The following default values 
have been assigned to the elements of CORFAT: CORFAT(K) *> 0.2, 0.40, 

0.60, 0.80, and 1.0 for K • 1 to 5, respectively. 

ITC0R A positive Integer denoting the Iteration number at which the 

potential function under-relaxation is terminated. For iteration 
numbers greater than ITC0R, the standard algorithm with no 
under-relaxation is applied. Generally, fastest convergence is 
obtained by using under-relaxation for all iterations. The default 
and recommended value for ITC0R is 10,000. 

VJT A real variable denoting the exhaust jet velocity normalized by the 

critical sonic speed. A default value of 0.9 is specified for VJT. 

JNPR An integer which controls the meridional station output frequency 

when the nacelle surface solution is printed. Setting JNPR**l prints 
every circumferential station. A default value of 3 is specified for 
JN’PR. 

XDUMP A one-dimensional integer array consisting of seven elements. Each 
element of KDUMP specifies whether or not a particular nacelle 
algorithm subroutine is to have debug output printed. Specifying 
RDUMP(I)»l (I«l,...,7) activates the debug output option for the 
corresponding subroutine. Specifying KDUMPd) *■ 0 causes the debug 
printout to be skipped. The eleaents of KDUMP activate the debug 
printout for the following subroutines/conputations and have the 
following default values. 


K.DUMP(I) Activates Debug Output Default 
for (SUBROUTINE) Value 

RDUMP(i) metric -alculations (METRIC) 0 
RDUMP ( 2 ) Initialization (NINIT) 0 
KDUMP(3) not presently used 0 
KDUMP(4) AF2 solution scheme (NS0LVS) 0 
KDUMP(5) physical density calculation (NR0) 0 
RDUMP(6) residual calculation (NRESID) 0 
RDUMP(7) modified density calculation (h'ROCO) 0 


tTSTRT An integer denoting the iteration number at which debug output io to be 
initiated. A default value of 1 is specified for 1TSTRT. 

8.3 PARAMETER STATEMENT SPECIFICATION 

The TnN flow solution program uses variable array dimensions. The 
respective array sizes are fixed by specification of certain parameters on 3 
PARAMETER statement. The PARAMETER statement for the TWN code is identical to 
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that used in tha LK grid interface code. Hence, the user should refer to 
Section 6.2 of thia chapter for a definition, of the dimension paraneters for 
the TWH code. 


8.4 FILE USAGE 


The following files are used by the TWN flow solution code. 


File No. 


Usage 


TAPE 5 
TAPE 6 
TAPES 9,10 


TAPE U 


TAPE 12 


TAPS 19 


TAPE 20 


TAPE 44 


input file 

printed output file 

input files from the GRGEN3 wing grid 

generation code containing the wing grid 

coordinates 

input file containing the potential field 
from a previous wing solution (for restart 
option) 

output file on which the current wing 
potential field is written (for restart 
option) 

input file from the LK grid interface code 
containing the grid overlap region 
parameters 

output file containing the wing surface 
pressure coefficient field 
input file from the LK grid interface code 
containing the scaled nacelle grid coordinates 


9. T.N PROGRAM OUTPUT INTERPRETATION 


The initial portion of the computer printout for the TWN program consists 
of the NAMELIST input data. Then the preliminary information listed below is 
output. 


RHOSTR 

RHOINF 

QINF 


NI 


NJ 


NK 

IHALF 


critical density normalized by stagnation density 
free-stream density normalized by stagnation density 
free-stream velocity magnitude normalized by critical 
sonic speed 

maximum number of wing grid points in the wraparound 
direction 

maximum number of wing grid points in the spanwise 
direction 

maximum number of wing grid points In the normal direction 
equals NI/2+1 


The wing algorithm output parameters are defined below. 


ITER 

ICMAX 

JCMAX 

KCH.A w 

CMAX 

IRMAX 


iteration number 

location of maximum correction (£-index) 
location of maximum correction (q-index) 
location of maximum correction (r-index) 
maximum correction 

location of maximum residual (^ -index) 
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JRMAX 

KBMAX 

RMAX 

NSP 

ALPHA 

GAM 

in 

XU 

ZU 

XL 

ZL 

XUC 

XLC 

CPU 

CPL 

RHOU 

RHOL 

MACHU 

MACKL 

I 

Y 

Y/TIP 

CL/C 

CD/C 

CM(RTC/4) 


location of maximum rcaidual (n-index) 
location of nasiaua rcaidual (£-index) 
maxima residual 
nuabar of suparconic point3 
acceleration parameter 

an array containing circulation values for four consecutive 
iterations 

grid point nunbar on wing surface 

x-coordinate for wing upper surface 

z-coordinate for wing upper surface 

x-coordinato for wing lower surface 

z-coordinate for wing lower surface 

x-coordinate for wing upper surface in X of chord 

x-coordinate for wing lower surface in X of chord 

pressure coefficient on wing upper surface 

pressure coefficient on wing lower surface 

density on wing upper surface normalized by stagnation 

density 

density on wing lower surface normalized by stagnation 
density 

Mach number on wing upper surface 
Mach number on wing lower surface 
wing span station number 
y-coordinate of wing span station 

y-coordinate of wing span station in terms of fraction of 
wing seaispan 

local wing section lift coefficient scaled by local chord 
local wing section wave drag coefficient scaled by local 
chord 

local wing section moment coefficient about the root 
quarter chord 


The nacelle algorithm output parameters are defined below. 


I 

J 

X 

Y 

Z 

DX 

RHO 

MACH 

CP 

MSP 

RMAX 

RAVG 

I RMAX 

JRMAX 

KRMAX 


£ -wraparound pc.nt index 

n-circuaferential point index 

x-coordinate of surface solution point 

y-coordinate of surface solution point* 

z-coordinate of surface solution point * 

axial (x) distance between surface solution point and 

hilite point for a given meridional plane 

density coefficient defined as the local density 

normalized by stagnation density 

Mach number 

pressure coefficient 

number of supersonic points 

maximum residual 

average residual 

^-location of maximum residual 

H-location of maximum residual 

^-location of maximum residual 


* (Roto - the y and z axes are along the wingspan for the wing and 
nacelle algorithms, respectively; the z and y axes are normal to 
the wing for the wing and nacelle algorithms, respectively) 
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SAMPLE CASES AND SUGGESTIONS FOR USAGE 


1. INTRODUCTION 

In this section, sample cases are presented to illustrate application of 
both the GRGEN3 and NGRIDA grid generation programs, the LK grid interface 
program, and the TUN flow analysis program. For each sample case, a 
discussion of the problem is given, the required input data are presented, and 
selected portions of the computer output are listed. The input parameter 
discussions follow the order in which the input parameters are presented in 
Section IV. At the end of this section, suggestions are provided for use of 
the computer programs. 

2. GRGEN3 WING GRID GENERATION PROGRAM SAMPLE CASE 

This sample case is presented to illustrate application of the GRGENJ 
wing grid generation program to determining the computational grid about a 
NACA 0012 unswept and untapered wing configuration. This case represents the 
case being considered when the default wing geometry contained within the 
GRGEN3 code is used. 

The data deck for this sample case is Illustrated in Figure 35. All 
input parameters in NAMELIST GRID3D retain their default values except for NJ, 
NK, NJT, NJC, DYTIP, DYNCL, and DYR00T. The number of spanwise and 
normal stations to be used in the grid are 30 and 22, respectively; hence, NJ 
* 30 and NK * 22 are entered in NAMELIST GRID3D. The spanwise indices 
denoting the first station off the wing tip and the nacelle location are 25 
and 13, respectively; hence, NJT » 25 and NJC » 13 are specified. The 
spanwise spacings near the wing tip and wing root are 0.4 and 0.2, 
respectively; hence, DYTIP ■ 0.4 and DYR00T » 0.2 are specified. The spanwise 
grid spacing near the nacelle symmetry plane is controlled by entering 
DYNCL-0.20. 

All input parameters in NAMELIST GRIDIN retain their default values 
except for the grid point distribution parameter BINN. BINN ■ 1.08 is entered 
which clusters points near the leading and trailing edges of the wing surface. 

Selected portions of the computer output for this sample case are 
presented in Figure 36. 

3. NGRIDA NACELLE GRID GENERATION PROGRAM SAMPLE CASE 

This sample case is presented to illustrate application of the NGRIDA 
nacelle grid generation program to determining the interface computational 
grid about the Lockhepd-Georgia GULACi axisymmetric nacelle/ inlet 
configuration. This case represents the case being considered when the 
default nacelle geometry contained within the NGRIDA code is used. 

The data deck for this sample case is illustrated in Figure 37. All 
input parameters in NAMELIST LIST1 retain their default values except for 
KMAX, RAD0UT , XLEFT, ALPHAO, ALPHAI, and ITAPE. Five radial stations will be 
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SGRID3D HJ-30, NK-22, NJT-25, NJC-13, DYTIP-0.4, 
DYNCL-0.20, DYROOT-O.2 SEND 
SGRIDIN BXNH-1.08 SEND 


Figure 35. Data deck for GRGEN3 sample case 



I 


GRGEN3, THilEE DIMENSIONAL MING GRID GENERATOR USING A TWO DIMENSIONAL GRID AT EACH SPAN STATION. 

WRITTEN BY TERRY HOLST (APPLIED COMPUTATIONAL AERODYNAMICS BRANCH ) > ALAN FERNQUIST, AND SCOTT D. 

*** NASA AHES RESEARCH CENTER, MOFFETT FIELD, CALIFORNIA 9403S *** 

THE MAXIMUM DIMENSIONS ARE 151 BY 30 BY 30. 

SGRID3D 


DYROOT 

s 

0,2000000000000000 , 

DYTIP 

a 

0,4000000000000000 , 

DYNCL 

o 

0,2000000000000000 , 

MI 

a 

101, 

NJ 

a 

30, 

HJC 

a 

13, 

NJT 

□ 

25, 

NK 

a 

22, 

SFACT 

a 

1.240000000000000 , 

SMOOTH 

3 

l.OOOGCOOOOOOOOOO , 

SWPANG 

a 

0,00 0000 OOOOOOOOOOE+QO, 

TIP 

a 

6.000000000000000 , 

TWANG 

S3 

O.OOOOOOOOOOOOOOOOE+OO, 

CEL 

CEND 

a 

3.000000000000000 


NOTE THAT THERE ARE 101* 30= 3030 POINTS IN EACH K-PLANE, AND 101* 30* 22*» 66660 TOTAL POINTO, 
INCIDENTALLY, 101* 30* 22*6= 399960 


INTERNALLY GENERATED Y DISTRIBUTION 


FOR TIP AT Y =6,0000, AND HJT * 25 


J Y ( J ) 


1 

2 

3 

4 

5 


0.000000 

0,200000 


0.411177 

0,647392 

0,911773 



C 


Figure 36. Output for GRCEN3 sample case. 
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27 

7,000000 

28 

7,656000 

29 

8,347200 

30 

9,176640 


sgridin 


AHGRID 

a 

1.0000000000000000E-02 

ALGRID 

a 

6, 0000000000000000 S"05 

BINH 

9 

1.000000000000000 

ERGRID 

a 

9,S999999999?39999E®04 

IIHCR 

a 

3# 

IQPT 

a 

X, 

JIHCR 

a 

6# 

KGRID 

a 

0 , 

HAXIT 

a 

200, 

HGRID 

a 

8, 

NOUT5 

s 

0, 

NOUT6 

0 

201, 

NSTATN 

s 

0, 

OHEG 

c 

2.000000000000000 

RADHAX 

a 

6,000000000000000 

TMAX 

8 

0,1200000000000000 

XC 

a 

Oo'iOOOOOOOOOOOOOOO 

XCN 

a 

O.OOOOOOOOOOOOOOOOE+OO 


6END 

NOTE THAT FOR THE HACAOQXX OPTION THE THICKNESS-TO-CHORD RATIO IS ACTUALLY 0.1109, 


Figure 36. Continued 
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*** FINAL AIRFOIL COORDINATE DISTRIBUTION , Y * 0,000000 *** 


I 

XB 

Y b 

S 2 

SX 

I 

XB 

YB 

S 2 

SX 

I 

XD 

YQ 

1 

1.00000 

0.00000 

0.00000 

0.00267 

35 

0.21523 

0,00000 

0,78477 

0.02540 

69 

0.2696$ 

0.00000 

2 

0.99467 

0.00000 

0.00533 

0.00567 

36 

0.19071 

0.00000 

0.00929 

0.02364 

70 

0.29501 

0.00000 

3 

0.90866 

0.00000 

0.01134 

0, 00639 

37 

0.16795 

0,00000 

0.03205 

0.02169 

71 

0.33003 

0,00000 

4 

0.90189 

0.00000 

0,01811 

0.00710 

38 

0,14694 

0.00000 

0,05306 

0.02015 

7 a 

0.3623$ 

0,00000 

5 

0.97430 

0.00000 

0,02570 

0,00806 

39 

0.12764 

0.00000 

0,07236 

0.01347 

73 

0.33376 

0.00000 

6 

0.96578 

0.00000 

0.03422 

0,00902 

40 

0.11000 

0,00000 

0.89000 

0.01635 

74 

0.43003 

0,00000 

7 

0.95625 

0.00000 

0,04375 

0.01000 

41 

0,09395 

0.00000 

0.90605 

0,01531 

75 

0.66406 

0,00000 

8 

0.94561 

0.00000 

0,05439 

0.01124 

42 

0,07939 

0,00000 

0,92061 

0.01266 

76 

0.50000 

0.00000 

9 

0.93376 

0,00000 

0,06624 

0,01250 

43 

0.06624 

0,00000 

0,93376 

0.01250 

77 

0,53514 

0.00000 

10 

0,92061 

0.00000 

0,07939 

0,01386 

44 

0.05439 

0,00000 

0,94561 

0.01124 

79 

0.56993 

0,00000 

11 

0,90603 

0,00000 

0,09395 

0,01531 

45 

0.04375 

0,00000 

0.95625 

0.01000 

79 

0,60624 

0,00000 

12 

0.89000 

0.00000 

0.11000 

0.01695 

46 

0,03422 

0.00000 

0,96578 

0,00902 

80 

0.63765 

0,00000 

13 

0.87236 

0.00000 

0,12764 

0,01847 

47 

0.02571 

0,00000 

0.97430 

0.00006 

01 

0*06397 

0,00000 

14 

0.85305 

0,00000 

0,14694 

0,02015 

40 

0.01009 

0,00000 

0.90109 

0,00718 

02 

0.70089 

0,00000 

15 

0.83205 

0.00000 

0.16795 

0,02189 

49 

0.01139 

0,00000 

0.90066 

0.00639 

03 

0,73054 

. 0.00000 

16 

0.00929 

0,00000 

0.19071 

0.02364 

50 

0.00523 

0,00000 

0,99467 

0.00567 

84 

0,75050 

0,00000 

17 

0.78477 

0.00000 

0.21523 

0,02540 

51 

0. 00000 

0,00000 

1.80000 

0.00533 

05 

0,73477 

0.00000 

18 

0,75850 

0,00000 

0.24150 

0.02711 

52 

0.00523 

0.00000 

1.00533 

0,00567 

06 

0.00929 

0.00000 

19 

0,73054 

0,00000 

0,26946 

0,02076 

53 

0.01139 

0,00000 

1.01134 

0.00639 

87 

0.03303 

0.00000 

20 

0.70099 

0.00000 

0,29901 

0,03029 

54 

0.01809 

0.00000 

1.01811 

0.00718 

00 

0.0530$ 

0,00000 

21 

0.56997 

0,00000 

0,33003 

0,03167 

55 

0.02571 

0,00000 

1.02570 

0,00306 

09 

0,07236 

0,000:0 

22 

0.63765 

0,00000 

0,36235 

0,03286 

56 

0.03422 

0.00000 

1.03422 

0,00902 

90 

0. Q 9000 

0.00000 

23 

0.60424 

0.00000 

0,39576 

0,03383 

57 

0.04375 

0,00000 

1.04375 

0.01008 

91 

0,90605 

0.00000 

24 

0,56998 

O.OOOOu 

0,43002 

0,03455 

50 

0.05439 

0,00000 

1.05439 

0.01124 

92 

0.92061 

0.00000 

25 

0.53514 

0,00000 

0.46486 

0,03499 

59 

0,06624 

0,00000 

1.96624 

0,01250 

93 

0.93375 

0.00000 

26 

0,50000 

0,00000 

0.50000 

0.03514 

60 

0.07939 

0,00000 

1,07939 

0.013 S 6 

94 

0.96561 

0.00000 

27 

0.46406 

0.00000 

0.53514 

0.03499 

61 

0,09395 

0,00000 

1.09395 

0.01531 

95 

0.95625 

0.00000 

28 

0.43002 

0.00000 

0.56990 

0.03455 

62 

0.11000 

0,00000 

1.11000 

0,01665 

96 

0.96570 

0,00000 

29 

0.39576 

0.00000 

0.60424 

0.03383 

63 

0.12764 

0.00000 

1,12764 

0. 01047 

97 

0,97430 

0,00000 

30 

0.36235 

0.00000 

0.63765 

0.03286 

64 

0.14694 

0.00000 

1,14694 

0.02015 

90 

0.93139 

0,00000 

31 

0.33003 

0.00000 

0.66997 

0.03167 

65 

0.16795 

0,00000 

1.16795 

0.02109 

99 

0,90066 

0.00000 

32 

0.29901 

0.00000 

0,70099 

0.03029 

66 

0.19071 

0,00000 

1.19071 

0,02364 

100 

0,99467 

0,00000 

33 

0,26946 

0.00000 

0.73054 

0.02376 

67 

0.21523 

0,00000 

1.21523 

0.02540 




34 

0.24150 

0,00000 

0.75050 

0.02711 

68 

0.24150 

0,00000 

1,24150 

0.02711 






Figuie 36. Continued. 
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$LIST1 KllAX-5, RAD0UT-16.0, 
ALPHAI-1.09, ITAPE-14 $END 
SLIST2 $E»D 
SHST3 SEND 
SLIST4 SEND 


XLEFT— 12.0, ALPHAO-1.09, 


Figure 37. Data deck for NGRIDA sample case. 
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uacd in the interface grid, hence KMAX“5 is specified in NAMELIST LIST1. 
Because an interface grid is being generated, the magnitudes of the outer 
boundary parsmetoro RADOUT and XLEFT are reduced from their default values to 

16.0 and -12.0, respectively. The ALPHAO and ALFHAI grid stretching function 
parameters are changed slightly from their default values to a value of 1.09 
each. Since this is an interface grid generation, the nacelle grid coordinate 
data oust be loaded onto disk file TAPE 14; hence, ITAPE “ 14 is specified. 

All input parameters entered in NAMELISTS LIST2, LIST3, and LIST4 retain 
their default values. 

Selected portions of the computer output for this sample case are 
presented in Figure 38. 

4. LK GRID INTERFACE PROGRAM SAMPLE CASE 

This sample case is presented to illustrate application of the LK grid 
interface program. The case considered interfaces the NACA 0012 wing grid and 
the GELACl nacelle grid generated by the sample cases of the GRGEN3 and NGRIDA 
programs, respectively. 

The data deck, for this sample case is illustrated in Figure 39. All 
input parameters in NAMELIST NACEL retain their default values except for XR. 
The distance, expressed in terms of the fraction of the local chord, between 
the nacelle hilite point 3nd the wing leading edge point is 1.0. Hence, XR » 

1.0 is specified in NAMELIST NACEL. 

Selected portions of the computer output for this sample case are 
presented in Figure 40. 

5. TWN FLOW SOLUTION PROGRAM SAMPLE CASE 

This sample case is presented to illustrate application of the TWN flow 
analysis program. This case considers the multicomponent flow field 
computation about the NACA0012-wing/GELACl-nacelle configuration. This case 
employs the grids and interface parameters generated by the GRGEN3, NGRIDA, 
and LK sample case executions. 

The data deck for this sample case Is illustrated in Figure 41. All 
input parameters in NAMELIST FL0WIN retain their default values except for 
ALPHAW, NDIF, RGAM, and JSKIP. For this case, the free-strean Mach number is 
retained at its default value of 0.75. The angle of attack considered in this 
case equals 1.0 degree, hence ALPHAW"1.0 is entered in NAMELIST FL0WIN. 
Upwinding In regions of supersonic flow will be applied only in the wraparound 
direction, hence NDIF=0 is entered. The wing circulation relaxation parameter 
RGAM is entered as 1.1. The spanwise output frequency for the wing solution 
is specified as JSKIP-12. A title card is entered after NAMELIST FL0UIN. 

All input parameters in NAMELIST FL0WC retain their default values except 
for CRATI0 and JNPR. For this case, the inlet capture ratio CRATI0 is 
specified as 0.7. The circumferential output frequency for the nacelle 
solution is now specified as JNPR"6. 


Selected portions of the computer output for this sample case are 



TmR EE- DIMENSIONAL NACELLE GRID GENERATION PROGRAM 


LISTING 

OF 

INPUT DATA 


9LIST1 

1MAX 

3 

67, 


JMAX 


13, 


KMAX 

3 

5, 


XLEFT 

3 

-12.00000 

9 

RADOUT 

o 

16.00000 

9 

DELTA 

8 

0.5000000 

9 

IOUTER 

3 

30, 


I INNER 

3 

30, 


ALPHA I 

3 

1.09000G 

9 

ALPHAO 

a 

1.090000 

9 

KQUTER 

s 

0, 


ALPHB1 

s 

1.00000U 

9 

ALPH30 

3 

1.000000 

9 

SCALE 

3 

0.1190476 

9 

ITAPE 

3 

14, 


KPLOT 

M 

1, 


KDEV 


1, 


KDUHP 

= 

1. 


IAA 

c 

12, 


IDD 

3 

57, 


KPRINT 

8 

1. 

24*0 

KTIME 

3 

3, 


NDEL 

3 

5, 


DSFACT 

Z 

1.000000 

9 

KRAMP 

8 

1. 


IRAHP 

3 

57 


SEND 

SLIST2 

NO 

3 

100, 


XOUT 

3 

0.0000000£+00, 


4.9999999E-03, 9, 9999998E-03 , 1.5000000E-02, 


2.0000000E-02, 
4. 5000002E-02 , 
7. GOGOOGOE-02 , 


2 , 5000000E-02 » 
5.00U0001E-02, 
7 • 9999&90E-O2 , 


2.9999999E-02, 

5.5000000E-Q2, 

9.0000004E-02, 


3.5000000E-02, 
5.99S9999E-02, 
0. 1000000 , 


3.9999999E-02, 
6 »4999990E-O2 , 
0.1100000 , 


Figure 38. Output for NOKIDA sample case. 
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GRID POINT COORDINATES FOR SELECTED MERIDIONAL PLANES 


la 

1 

J* 

1 

K= 

1 

Xa 

0.619048E+01 

Ya 

0. 190476E+01 

Za 

O.OOOOOOt+OO 

I a 

2 

J3 

1 

K* 

1 

Xo 

0.581150E+01 

ys 

0, 190476E+01 

Z = 

O.OOOOOOE+OO 

1 = 

3 

J3 

1 

Ka 

1 

Xa 

0 . 544038E+0 1 

Ya 

0. 190476E+01 

Z= 

0.000000^+00 

1 = 

4 

J3 

1 

Ka 

1 

Xr 

0.SO9019E+01 

y= 

0.1S0476E+01 

2a 

0.000000^+00 

1 = 

5 

J* 

1 

K = 

1 

Xa 

0.47667OE+01 

y= 

0. 19U476E+01 

Za 

0 , 000000&+00 

1 = 

6 

Ja 

1 

Ks 

1 

Xa 

0.446923E+01 

Y= 

0.190476E+01 

Za 

0 . OOOQOOE+OO 

1 = 

7 

J = 

1 

Ka 

1 

X a 

0.419443E+01 

ya 

0. 190476E+01 

Za 

o . ooooool+oq 

1 = 

e 

Ja 

1 

K = 

1 

Xa 

0.393825E+01 

y= 

0. 190476 E + 01 

Za 

o.ogooooe+oo 

1 = 

9 

Ja 

1 

Ka 

1 

Xs 

0 , 36947 3E+0 1 

Ya 

0.190476E+01 

2 = 

o.oooooos+oo 

1 = 

10 

J = 

1 

Ka 

1 

Xa 

0.345491E+G1 

Ya 

0, 15047oE+01 

z= 

G.OOOOOQE+OO 

1 = 

11 

J = 

1 

Ka 

1 

Xa 

0.32G696E+01 

Ya 

0, 190476E+01 

Za 

O.OUOOOOE+OO 

1 = 

12 

Ja 

1 

K= 

1 

X = 

0.294156E+01 

Ya 

0,19047 6E+0 1 

Z= 

O.OOOOOOE+OO 

1 = 

13 

Ja 

1 

K = 

1 

Xa 

0.266812E+01 

Ya 

0, 19G476E+01 

Za 

O.OOOOOOK+OO 

1 = 

14 

Ja 

1 

K= 

1 

Xa 

0 . 241 706E+0 1 

Ya 

0.190476E+01 

Za 

O.OOOOOOL+QO 

1 = 

15 

Ja 

1 

Ka 

1 

Xa 

0.22077 1E+01 

Ya 

0. 19O476E+01 

Za 

0. 000000^+00 

1 = 

16 

Ja 

1 

K = 

1 

Xa 

0.203231E+01 

Ye 

0. 19047bE+01 

Za 

o.oooooo^+oo 

1 = 

17 

Ja 

1 

K = 

1 

Xa 

0.1H7844E+01 

Ye 

0.190476E+01 

Za 

0.000000E«00 

1 = 

IB 

J = 

1 

Ka 

1 

Y.r. 

0. 175016E+01 

Ya 

0. 190476E+01 

Z= 

0.O00000E+00 

1 = 

19 

Ja 

1 

K = 

1 

Xa 

0 . 163354 Et 01 

Ya 

0 . 19047 oEtG 1 

Za 

O.OOOOGO^+OO 

1 = 

20 

J a 

1 

Ka 

1 

Xa 

0. 152214E+01 

Ya 

0. 19u47bE+01 

Z = 

0 . OOOOOOK+OO 

1 = 

21 

Ja 

1 

Ka 

1 

Xa 

0.1413 0 3E + 01 

Ya 

0. 190476E+01 

Z = 

0 , 000000^+00 

1 = 

22 

Ja 

1 

Ka 

1 

Xa 

0. 1303S4E+01 

Ya 

0. 190476E+01 

Za 

0 , 000000^+00 

1 = 

23 

Ja 

1 

Ka 

1 

Xa 

0, 119096E+01 

Ya 

0, 190476Et01 

Za 

0 . 000000^+00 

1 = 

24 

Ja 

1 

Ka 

1 

X = 

0. 1O7226E+01 

Ya 

0. 190476E+01 

Z= 

0.000000ti + 00 

Ic 

25 

J = 

1 

Ka 

1 

Xa 

0 , 943739E+00 

Ya 

0. 140476E+01 

Za 

0 . 000000^+00 

1 = 

26 

Ja 

1 

Ka 

1 

Xa 

0 . B00539E+00 

Ya 

0. 190476E+G1 

Z= 

o.nooooo^+oo 

1 = 

27 

Ja 

1 

Ka 

1 

X» 

0.635786E+00 

Ya 

0.1 90476E+0 1 

Za 

0.000000t+00 

1 = 

28 

Ja 

1 

Ka 

1 

X s 

0 . 439 1 WbE+00 

Y = 

0.19047ot>01 

2 = 

O.OGOOOOt'+OO 

Ic 

29 

J = 

1 

Ka 

1 

Xa 

0.251400E+00 

Ya 

0. 18U206E+01 

Zs 

O.OOOOOOL+OO 

1= 

30 

Ja 

1 

Ka 

1 

Xa 

0.1091&0E+C0 

Ya 

0.179077 E+0 1 

Za 

0.0o0000t+00 






Figure 38. Continued 
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$NACEL XR-1.0 $END 


Figure 39, Data deck for LK sample case. 
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ff. i > ■ « / ■— i *» : ,' 5 ^ 


ssEgas s^vi ^agre s sjgTi 32sasj'= 


vO 

v_n 


5NACEL 


XR 

a 1.000000 

ZT 

a 0.4000000 

SCALE 

a 0,2200000 

SEND 



XR* *1» 00000 XLo 0,04762 ZT« “0,40000 Zfic -0,87143 
XWORa -1,00000 YWORo 3,00000 ZWORe -0,63571 
NJCal3 Jl=ll J2 b15 

NUMBER OF POINTS ON ZIGZAG BOUNDARY N = 44 

Kl= 9 K2= 17 

NUMBER OF POINTS ON ZIGZAG SURFACE NS e 94 

NO. OF EXHAUST POINTS NN r 5 
XJls 14 KJ2«= 16 

Figure 40. Output for LK sample ci.se. 


O O 
12 




SFLOWIN ALPHAW-1 .0, NDIF-0, RGAM-1.1, JSKIP-12 SEND 
NACA0012/GELAC1 
SFLOWC CRATIO-O.7, JNPR-6 SEND 


Figure 41. Data deck for TWN sample case. 
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presented in Figure 42. 

6. SUGGESTIONS FOR USAGE 
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The CRGEH3 and NCR IDA grid generation program, the LK grid interface 
program, and the TON flow analysis program havo been executed on the 
VAX-1 1/780 and CRAY-l conputcra. All of these programs employ PARAMETER 
statements which fix the respective array dimensions at the tine of 
compilation. This provides an effective means of changing program core 
storage requirements. The VAX version of each program employs an INCLUDE 
statement which inserts the respective PARAMETER statement at the appropriate 
locations within the code. The GRGEN3 wing grid generation progran employs 
the 'INCLUDE 'GRGEN3.DIM" statement. The NGRIDA nacelle grid generation 
program employs the 'INCLUDE 'NGRID.DIM" statement. The LK grid interface 
program employs the 'INCLUDE 'TWING.DIM'' and 'INCLUDE 'NACELLE. DIM" 
statements. The TWN flow analysis program employs the same INCLUDE statements 
as does the LK program. If the codes are executed on a VAX, then the 
GRGEN3 .DIM, NGRID.DIM, TWING.DIM, and NACELLE. DIM files must reside in the 
user's directory. For execution on the CRAY-1 or CDC computers, the PARAMETER 
statement is edited Into the respective code by use of a COMDECK. For use on 
other computer systems, the PARAMETER statements will have to be inserted in 
the respective programs at the locations noted by the INCLUDE statements. 
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ft 


8FL0WIN 




ALPHAH 

a 

4.000000 , 


ALPHAL 

X 

0,4000000 , 


ALPhAW 

a 

1.000000 , 


ARAT 

a 

0 , OOGOOOOE+OO, 


BXI 

a 

0,1000000 , 


CON 

a 

1,200000 , 


ERROR 

a 

1 . OOOOOOOE-03 , 


FMACH 

a 

0.7500000 , 


IfITOUX 

a 

0, 


IOUT 

a 

4, 


JCPEHD 

a 

0, 

o o 

JSKIP 


12, 

JO 

KKK 

a 

0, 

~a £ 

o 5 

MAXIT 

a 

40, 

o > 

M 

a 

8, 

■a r 

MCP 

= 

o, 

<o -o 

C V 

MORIN 

s 

0, 

> o 

MOROUT 

= 

0, 

C r * - 

NDIF 

3 

0, 

- J c a 

OMEGA 

S 

1,600000 , 


RGAM 

3 

1,100000 , 


TIP 

a 

0 , OOOOOOOE + OO , 


MODE 

a 

1, 


KCX 

SEND 

- 

4 



Figure 42. Output for TWN sample case. 



3D TRANSONIC FULL POTENTIAL WING-NACELLF. PROGRAM WITH MAXIMUM DIMENSIONS 121 BY 30 BY 25, 

NACA001 2/GEL AC 1 

Nla 101 NJa 30 NKs 22 NJTo 25 


FHACH 

ALPHAW 

ALPHAH 

ALPHAL 

RHOSTR 

GAMMA 

0.750 

1.000 

4.000 

0,400 

0,634 

1.400 

BXI 

OMEGA 

RHOINF 

OINF 

NJT 


0.100 

1,800 

0.766 

0.77y 

25 


N I 

NJ 

NK 

IHALF 



101 

30 

22 

51 




Figure 42. Continued 
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TRANSONIC HACELLE/INLET FLOW VARIABLES 


LISTING 

OF 

INPUT DATA 


8FL0WC 

XMCF 

3 

0.3516000 

9 

CRATIO 

8 

0,7000000 

9 

KINIT 

3 

1, 


UMAX 

3 

50, 


CHIT 

3 

1.0000000E-03 

9 

XMULM 

3 

0,7000000 

9 

KALPHA 

3 

1 » 


ALPHAL 

3 

1.000000 

9 

ALPHAH 

3 

20.00000 

9 

MO 

3 

Q. 


OMEGA 

8 

1,000000 

9 

NKVIS 

3 

0, 


CFACT 

3 

1.200000 

9 

BX IE 

3 

0.1000000 

9 

KSMTH 

3 

4, 


NOEL 

3 

3, 


KTIME 

3 

3. 


VJT 

3 

0.9000000 

9 

ITCOR 

3 

10000, 


CORFAT 

3 

0.2000000 

, 0 f 4000000 

1.000000 

, 5*0.0000000£+00, 

KDUMP 

C 

20*0, 


ITSTRT 

3 

1. 


JNPR 

cend 

S 

6 



0.6000000 


0.8000000 


Figure 42 


Continued 
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SPAN STATION ¥(13) «* 3.000 IS 50.00 PERCENT UF SEMISPAN 


IU 

XU 

ZU 

XL 

ZL 

fcXUC 

fcXLC 

CPU CPL 

RHOU 

khol 

HACHU 

MACtiL 

50 

0.000 

0,003 

0.000 

-0.003 

0,0 

0,0 

1,04533 1.14357 

0.9799 

0.9990 

0,2019 

0,0444 

49 

0.002 

0.008 

0.002 

-0.008 

0.2 

0.2 

0.69369 0.96144 

0.9103 

0.9635 

0,4370 

0.2730 

48 

0.006 

0.013 

0,006 

-0,013 

0.6 

0.6 

0.26737 0,65560 

0.8270 

0,9026 

0,6203 

0.4574 

47 

0.012 

0,018 

0.012 

-0.018 

1.2 

1.2 

-0,02445 0.38247 

0,7608 

0.8468 

0,7602 

0,5865 

46 

0.019 

0.023 

0.019 

-0.023 

1.9 

1.9 

-0.22963 0,18044 

0.7159 

0.8045 

0,0457 

0,5741 

45 

0.027 

0.027 

0,027 

-0.027 

2.7 

2.7 

•0.36646 0.03341 

0.6853 

0.7732 

0,9032 

0,7360 

44 

0.036 

0.031 

0.036 

-0.031 

3.6 

3.6 

-0.46499 -0,07799 

0.6630 

0.7492 

0,9452 

0.7025 

43 

0.047 

0,034 

0.047 

—0.034 

4.7 

4.7 

-0,54266 -0.16625 

C.64bl 

0.7299 

0.9760 

0,0192 

42 

0, 0a9 

0.038 

0.059 

-0,038 

5.9 

5.9 

-0.61500 -0.23e57 

0.6284 

0.7139 

1,0105 

0.3494 

41 

0.073 

0.041 

0.073 

-0,041 

7.3 

7.3 

-0.66221 -0.29675 

0.6173 

0,7 005 

1,0315 

0,0746 

40 

0.088 

0.044 

0.086 

-0.044 

8,6 

8.8 

-0.70372 -0.34682 

0.6076 

0.6893 

1,0502 

0,0957 

39 

0.105 

0.047 

0.105 

-0.047 

10.5 

10,5 

-0,74127 -O.39030 

0.5987 

0.6800 

1.0672 

0.9133 

38 

0.124 

0.050 

0.124 

-0,050 

12.4 

12.4 

-0.77581 -0.42347 

0.5904 

0.6725 

1.0031 

0.9274 

37 

0.144 

0.052 

0.144 

-0.052 

14,4 

14.4 

-0.80010 -0,44924 

0.5846 

0.6666 

1.0943 

0.9384 

36 

0 . 1 06 

0,055 

0,166 

-0.055 

16,6 

16,6 

-0,81403 -0.46697 

0.5813 

0,6621 

1.1003 

0.9469 

35 

0.190 

0.056 

0.190 

-0.056 

19.0 

19,0 

-0.81722 -0.48340 

G.5805 

0.6500 

1,1023 

0.9531 

34 

0.216 

0.056 

0.216 

-0.058 

21,6 

21,6 

-0,80697 -0.48966 

0,5630 

0.6574 

1,0975 

0,9558 

33 

0.244 

0.Q59 

0.244 

-0,059 

24.4 

24.4 

-0,77543 -0,46279 

0.5905 

0,6539 

1.0029 

0,9528 

32 

0.273 

0,059 

0,273 

-0.059 

27.3 

27,3 

-0,70482 -0.46268 

0.6073 

0.6635 

1,0507 

0.9442 

31 

0.304 

0.059 

0.304 

-0.059 

30.4 

36.4 

-0,56330 -0.43365 

0.6404 

0.6701 

0,9070 

0.9319 

30 

0,336 

0,059 

0.336 

-0.059 

33.6 

33.6 

-0,50872 -0.40035 

0.6530 

0.6777 

0,9640 

0,9176 
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Figure 42. Continued. 
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29 0.370 0,033 0,370 -0,058 

20 0.404 0.0S7 0.404 -0.057 

27 0,440 0.056 0.440 -0,056 

26 0.475 0,054 0,475 -0.054 

25 0.511 0.051 0.511 -0.051 

24 0.547 0.049 0.547 -0,049 

23 0.532 0.046 0.532 -0,046 

22 0,616 0,043 0,616 -0,043 

21 0.649 0.040 0,649 -0,040 

20 0,602 0.037 0.602 -0.037 

19 0.712 0,034 0.712 -0,034 

10 0,741 0,032 0.741 -0,032 

17 0,7t>9 0.029 0,769 -0,029 

16 0.7y5 0.026 0,795 -0.026 

15 0.019 0.023 0.U19 -0.023 

14 0.041 0.021 0.841 -0.021 

13 0,861 0,018 0.861 -0,018 

12 0. 8 oO 0,016 0,880 -0.016 

11 0,897 0.014 0.897 -0.014 

10 0.912 0.012 0,912 -0,0:2 

9 0,926 0,010 0,926 -0.010 

0 0.939 0,000 0,939 -0.000 

7 0.950 0.007 0.950 -0.007 

6 0.961 0.005 0.961 -0.005 

5 0.970 0.004 0.970 -0.004 

4 0.978 0.003 0,978 -0.003 

3 0.985 0.002 0.985 -0,002 

2 0.992 0.001 0.992 -0.001 

1 0,997 0.000 0.997 0.000 

*** LIFT, DRAG 
CL a 0 , 144Q6BE+00 


37.0 37.0 -0.48208 -0.365oy 

40.4 40.4 -0.43954 -0.32972 

44.0 44.0 -0.39680 -0,29516 

47.5 47,5 -0.35581 -0.2616b 

51.1 51.1 -0.31670 -0,23022 
54.7 54,7 -0.27996 -0,20025 

58.2 58.2 -0,24523 -0.17102 

61.6 61,6 -0,21249 -0,14497 

64.9 64,9 -0.18155 -0.11942 

68.2 68,2 -0,15207 -0,09499 

71.2 71,2 -0,12406 -0,07158 

74.1 74.1 -0.097C5 -0,04890 

76.9 76.9 -0.07100 -0.02673 

79.5 79,5 -0.04553 -0.00493 

01.9 81.9 -0.02050 0.01603 

84.1 84.1 0.00425 0.03845 

86.1 86,1 0.02864 0,06039 

8e,0 88.0 0,05377 0,00265 

89.7 89,7 0,07698 0.10559 

91.2 91.2 0.10479 0.12936 

92.6 92.6 0.13173 0,15440 

93.9 93.9 0.16001 0.18061 

95.0 95,0 0.16992 0,20921 

96.1 96,1 0,22126 0.23938 

97,0 97,0 0.25279 0.26969 

97.8 97,8 0.27705 0.29448 

96.5 96.5 0.30290 0.31932 

99.2 99.2 0,32817 0,34423 

99.7 99.7 0.35343 0,36919 

AND MOMENT COEFFICIENTS FOR 
CO « 0 , 138965E-02 CM 


0.6591 0,6857 0.9523 0.9026 
0.6668 0.6936 0.9343 0.CC77 
0.6705 0.7013 0,9161 0.0731 
0,6677 0,7007 0,0907 0,0592 
0.6965 0.7150 0,0022 0.0459 
0.7047 0,7224 0,0560 0.0334 
0.7124 0.7207 0.0522 0.0216 
0,7197 0,7345 0.0305 0,0104 
0,7265 0,7401 0,e255 0.7297 
0.7330 0.7455 0,0133 0.7Q9S 
0,7391 0,7506 0,0017 0.779S 
0,7450 0,7555 0,7906 0,7704 
0,7507 0,7603 0,7796 0.7611 
0,7562 0,7650 0,7690 0.7S21 
0.7616 0,7697 0.7505 0.7430 
0.7670 0,7743 0.74Q2 0,7339 
0.7722 0,7790 0,7300 0,7243 
0.7776 0.7030 0,7275 0,7154 
0.7830 0.7037 0,7170 0,7050 
0,7865 0.7937 0,7061 0,6350 
0,7942 0,7990 0,6940 0.6Q52 
0.8002 0.0046 0,6320 0.6740 
0.8065 0, 6 106 0,6701 0,6619 
0.6131 0,bl69 0.6567 0,6490 
0.0190 0.8233 0,6432 0,6359 
0.6250 0,8265 0,6324 0,6252 
G ,o3G3 0.8337 0,6215 0.6143 
0.6355 0.8360 0,6105 0,6034 
0,6408 0.0440 0.5994 0.5924 

Y= 3,000 *** 

0.206606L-02 


TOTAL NUMBER OF CYCLES a 5 

Figure 42. Continued. 
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SECTION VI 

CONCLUSIONS AMD RECOMMENDATIONS 

An efficient grid interfacing zonal algorithm has been developed to 
compute the three-dlccnolonal transonic flow field about uing/nacelle 
multicomponent configurations. In the present approach, a component-adaptive 
grid embedding schema has been employed in which the global computational grid 
is composed of a series of overlapped component grids, where each component 
grid Is optimized for a particular geometry such as the wing or nacelle. The 
AF2 approximate factorization algorithm is used to determine the 
full-potential equation solution on each component grid with trivariate 
interpolation being used to transfer property information between the grids. 
The optimized grid system produces the desired balance between convergence 
speed and accuracy of the flow field solution. Numerical results using the 
overlapped grid scheme indicate that the present algorithm promises to be 
effective in computing the flow field about multicomponent configurations. 

To improve the accuracy of the present algorithm, it is suggested that a 
viscous correction capability be added to the inviscid flow prediction 
algorithm. This can be accomplished by either adding a boundary layer 
displacement thickness distribution onto the original body contours, or by 
using surface transpiration velocity boundary conditions in which the 
effective transpiration velocity is calculated from the computed boundary 
layer growth. Although the boundary layer flows on both the wing and nacelle 
are three-dimensional, using a two-dimensional boundary layer analysis applied 
in a stripwlse sense at span stations on the wing and at meridional stations 
on the nacelle may prove to be adequate for effecting a viscous correct ion to 
the inviscid algorithm. 

In order to correctly predict the uing lower surface pressure 
distribution near the wing leading edge, it will be necessary to complete the 
pylon modeling. To model a pylon of realistic thickness, it will probably be 
necessary to introduce a third component C-type grid around the pylon or to at 
least adjust the wing and nacelle component grids to have grid lines aligned 
with the pylon leading and trailing edges. 

It also appears that using an H-type wing grid may prove beneficial in 
attempting to interface the wing and nacelle flow solution algorithms. In an 
H-type wing grid, the mesh has greater uniformity and the mesh cells are more 
Cartesian-like underneath the wing surface in the vicinity of the nacelle. 

This should aid in interpolating flow properties from one component grid to 
another. 
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